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ABSTRACT 


A  practical  design  method  for  highly-loaded  blades  in  an  isolated 
cascade,ls  presented  in  this  thesis.^>The  flow  is  assumed  to  be 
incon^ressible  and  inviscid.  The  upstream  inlet  flow  condition  is  taken 
to  be  uniform.  The  £resen goals  of  this  research  are  to  provide  a 
practical  numerical  code  for  the  design  problem,  and  a  non-linear  theory 
which  can  be  easily  expanded  to  three-dimensions.  The  theory  is  based  in 
part  on  the  Clebsh  formulation.  The  blade  profile  is  determined 
iteratively  through  the  blade  boundary  conditions  using  a  "smoothing" 
technique.  A  practical  numerical  code  is  presented  for  the  design  problem 
using  "partial  smoothing".  The  program  gives  very  fast  convergence 
solutions  with  satisfactory  accuracy  for  practical  solidity  range. 
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Chapter  1  :  Introduction 

This  thesis  presents  a  technique  for  solving  the  design  problem  (the 
inverse  problem)  for  highly  loaded  blades  in  an  isolated  cascade.  In  the 
present  study,  the  flow  is  treated  as  if  it  were  incompressible  and  inviscid,  and 
the  upstream  inlet  flow  condition  is  assumed  to  be  uniform. 

Historically,  there  have  been  several  approaches  to  the  design  problem. 
In  one  approach,  for  example,  the  velocity  distributions  on  both  surfaces  of  the 
blade  are  specified,  and  the  resulting  blade  shape  is  calculated.  The  advantage 
of  this  technique  is  that  the  designer  can  prescribe  surface  pressure 
distributions  which  minimize  the  chance  of  flow  separation.  However,  the 
resulting  blade  geometry  is  not  guaranteed  to  have  a  realistic  configuration: 
the  blade  may  be  wavy,  or  even  have  negative  thickness. 

In  a  second  approach,  the  velocity  distribution  on  one  surface  and  the 
profile  thickness  distribution  are  specified,  and  the  resulting  blade  shape  is 
calculated.  This  formulation  does  not  guarantee  obtaining  an  acceptable 
pressure  distribution  on  the  second  blade  surface.  Moreover,  the  resulting 
loading  distribution  may  not  be  structurally  favorable:  for  example,  the 
loading  may  be  maximum  where  the  blade  is  thinnest. 

In  a  third  approach,  the  loading  and  thickness  (or  ’blockage’)  distributions 
are  specified,  and  the  resulting  blade  shape  is  calculated.  This  formulation 
does  not  guarantee  giving  acceptable  pressure  distributions  on  either  of  the 
blade  surfaces,  but  the  resulting  blade  shape  can  be  guaranteed  to  be  at  least 
structurally  sound. 
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In  applying  classical  aerodynamics  methods  to  these  problems,  the 
presence  of  the  blades  and  their  effects  on  the  flow  can  be  modeled  by 
distributing  singularities  (vortices,  sources  and  sinks)  on  the  blade  camber  or 
blade  surfaces.  Lewis  [lj  carried  out  the  design  problem  using  the  first  two 
approaches  mentioned  above.  His  analysis  is  based  on  Martensen’s  method  [2]: 
vortices  are  distributed  on  each  of  the  blade  surfaces,  and  the  induced  flow 
field  in  the  cascade  plane  is  calculated  using  the  Biot-Savart  law. 

Kashiwabara  [3j  carried  out  the  design  problem  using  the  first  of  the 
above  approaches  by  arranging  vortices,  sources  and  sinks  along  the  blade 
camber.  This  theory  also  attempts  to  take  into  account  some  three- 
dimensional  effects,  and  can  be  used  for  designing  blades  in  axial,  mixed  and 
radial  turbomachines. 

These  design  theories,  based  in  part  on  the  Biot-Savart  law,  are  simple 
for  cascade  calculation  only.  Until  recently  (4),  only  linearized  theories  have 
been  developed  to  design  three-dimensional  blades  (5],  (6j. 

One  goal  of  the  present  project  is  to  provide  a  practical  numerical  code 
for  the  design  problem  which  can  give  good  accuracy  and  fast  convergence 
solutions.  An  equally  important  goal  is  to  provide  a  non-linear  theory  which 
can  be  practically  expanded  to  three-dimensions.  The  present  theory  is  based 
in  part  on  the  Clebsh  formulation  which  has  been  successfully  investigated  by 
Tan,  Wang,  Hawthorne  and  McCune  [4]  in  their  study  of  a  three-dimensionel 
design  method  for  highly- loaded  but  infinitely  thin  blades. 

In  the  present  study,  the  blade  profile  is  determined  iteratively  through 
the  blade  boundary  conditions  using  a  "smoothing”  technique:  the  velocity 
potential  is  expressed  in  a  series  of  "smoothing”  functions  developed  by 


McCune  [Appendix  C].  The  "smoothing"  technique  represents  an  asymptotic 
expansion  (in  the  blade  spacing)  of  the  Green’s  function  for  blade  rows,  and 
can  be  applied  in  both  2D  and  3D. 

In  the  first  example,  the  blades  are  assumed  to  be  infinitely  thin  (the 
’Zero-Thickness’  problem)  to  show  the  power  of  the  "smoothing”  technique: 
the  blade  shape  is  solved  iteratively  through  a  set  of  algebriac  equations.  The 
results  compare  very  well  with  the  "exact”  method  [4].  Then,  a  similar 
approach  is  used  to  solve  the  inverse  problem  (in  the  third  formulation, 
outlined  above)  for  high  swirl  blades  having  prescribed  finite  thickness  (the 
’Finite-Thickness’  problem).  A  set  of  numerical  examples  are  presented,  which 
in  part  use  "partial  smoothing”,  i.e.  a  practical  truncation  of  the  smoothing 
series.  Blades  shapes  with  prescribed  thickness  (or  blockage)  distributions  are 
obtained  for  prescribed  swirl  schedules  or  loading  distributions.  To  confirm 
the  results  a  ’direct’  method  ,  developed  by  McFarland  [10],  is  used  to 
compute  both  the  circulation  and  the  pressure  coefficients  on  the  blade 
obtained  from  the  indirect  method.  The  results  show  that  the  blade  does  in 
fact  produce  the  desired  circulation  and  loading  distribution.  Moreover,  the 
pressure  distributions  on  each  surface  agree  well,  at  least  away  from  the 
leading  and  trailing  edges. 

Finally,  a  practical  design  procedure  is  presented  which  allows  for  rapid 
exploration  of  various  blade  loading  and  thickness  distributions,  so  as  to  be 
able  to  select  favorable  individual-surface  pressure  distributions  as  well. 


The  cascade  geometry  is  shown  id  the  above  figure.  The  flow  direction 
is  from  left  to  right.  All  lengths  are  non-dimensionalized  to  the  axial  chord 
length. 


The  blade  camber  lines  are  located  at: 


/$  =  spacing  between  blade  cambers 


f 


=  location  of  the  camber  line  n  =  0 


If  we  define  o(  as 

«  =  J  -  -foO 


then  blade  cambers  lie  on  surfaces  of  (X  =  rt/A  n  z  O,  i  I ,  i  2L  ... 

0  •  The  blockage  distribution 


Consider  the  blade  located  at  yj 
TOO  is  defined  such  that: 


-  the  blade  upper  surface  is  located  at  D  =  fOO  +  t(x) 

-  the  blade  lower  surface  is  located  at  j  ^  f(*)  -  T(x) 

Physically,  2T0C)  can  be  interpreted  as  the  axial  blockage 
distribution.  We  note  that  T(X)  is  not  the  blade  thickness  as  defined  in 

classical  aerodynamics,  where  the  thickness  is  defined  as  the  perpendicular 

distance  from  the  blade  camber  to  the  blade  surface. 


The  far  upstream  flow  condition  is  assumed  to  be  uniform.  All  velocities 
are  non-dimensionalized  to  the  upstream  inlet  x-component  velocity.  The  inlet 
flow  angle  is  denoted  by  Of,  .  Likewise,  the  outlet  flow  angle  is  denoted  by 


2.2  Fluid  mechanics  background 


As  mentioned  in  the  introduction,  we  have  chosen  to  model  the  presence 
of  the  blade  shape  by  distributing  singularities  on  the  blade  camber.  In  this 
section,  we  will  derive  the  governing  equation  for  the  velocity  field  in  the 
'  Zero-Thickness”  problem.  In  chapter  4,  we  will  extend  this  theory  to  the 
"Finite-Thickness”  problem. 


In  the  case  of  infinitely  thin  blades,  the  presence  of  each  blade  is 

modeled  by  arranging  vortices  along  the  blade  camber.  In  the  design  problem 

—  / 

the  swirl  schedule  V«(x)  is  given.  It  can  be  shown  that  the  swirl  schedule 
is  proportional  to  the  pressure  difference  (i.e.  the  loading)  accross  the  blade 
[Appendix  A], 

The  velocity  is  divided  into  two  parts:  an  average  velocity  y(*)  . 

rJ  . 

and  a  perturbed  velocity  >  *-e 

yc*o)  =  yoo  +  5(x,  y)  ,2.2-,, 

where  the  average  velocity  is  here  taken  to  be  the  pitch  average  defined  by: 

vc*)  =  —  j  vcm) 

and  thus,  by  the  definition  of  equation  (2.2-1),  we  must  have 

r ^ 

/  srcx,^ 

f 

The  above  equation  suggests  that  if  we  are  to  represent  <0*  as  a  series 


of  smooth  functions,  perhaps  we  should  choose  functions  which  possess  this 
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zero  pitch  average  property.  In  chapter  3,  we  will  show  that  the  "smoothing” 
functions  used  to  represent  &  possess  this  property. 

Consider  the  vorticity  field.  Given  the  assumption  of  incompressible, 
inviscid  and  uniform  inlet  condition,  then  by  Kelvin’s  theorem,  the  vortices 
must  lie  on  the  blade  camber,  and  the  flow  must  be  irrotational  everywhere 

else. 


It  can  be  shown  that  the  vorticity  is  related  to  the  blade  surface 
and  the  y-component  of  the  pitch  average  velocity  V  (*)  by  [Appendix 
B): 


-5-  (*/!?)  *  4  x  v&J 

where  Sp(?0  is  the  "periodic  delta"  function  [Appendix  B) 
Further,  G-(x)  =  -V^(X) 

The  pitch  average  vorticity  -XL.  is  defined  as: 


-£-60  E 


_  i 


C X/D-)  <^3 


=  W  x  VCr* 

The  vorticity  -XL.  is  defined  as: 

/V 

-a-  5  -5-. oo  -  4 

and  thus,  the  perturbed  vorticity  -XL  is 

A-  Oei'j)  =  -  ij  x  VG- 

-  V5(<0  X  VG 

where  S  (oc)  is  the  "sawtooth"  function  [Appendix  BJ. 


(2.2-2) 


(2.2-3) 


(2.2-4) 


We  are  now  in  the  position  to  write  down  the  governing  equation  for  the 
velocity  field  V  .  From  equation  (2.2-4),  we  can  write 

aJ 

Vx  £  -  VS  X  V(r  (2  2-S) 

The  Clebsh  formulation  says  that,  to  satisfy  (2.2-5),  write  the  perturbed 


velocity  as  as  the  sum  of  a  potential  part  and  a  rotational  part: 

£(*0)  =  V^<V,S)  +  &(x)VG-(x)  (2'2'6) 

rJ 

Note  that  the  curl  of  as  gives  _8_  back. 

A«  \l 

If  the  flow  is  incompressible,  then  V«  \1  —  0  =  V.  V  =  V*  IT .  Thus, 

V  V  /V 

from  equation  (2.2-6), 


v'-fic*,,)  =  -  v.  (svfr) 

=  -  t7S.  V&  -  S  Y% G- 


(2.2-7) 


This  gives  a  Poisson  equation  to  be  solved  for  <v>  with  appropriate 
boundary  conditions  (e.g.  on  the  blades). 


Finally,  the  velocity  field  in  the  cascade  region  is  of  the  form: 

V>o')  =  V(x)  +  fv$Z>(x,ij)  +  scx'iv&oo]  (2.2-8) 

In  the  design  problem,  V  is  given  and  we  seek  the  solution  for 

AS 

V  which  satisfies  all  necessary  boundary  conditions.  In  chapter  3,  we 

A/ 

present  a  design  method  for  infinitely  thin  blades  using  the  above  theory.  In 
chapter  4,  we  present  a  design  method  for  blades  having  finite  thickness  based 
in  part  on  the  above  approach.  In  this  case,  the  gap  average  velocity  is  given 
instead  of  the  pitch  average  velocity,  and  the  above  theory  needs  to  be 


Chapter  3  :  The  ’Zero-Thickness*  problem 


The  "Zero-Thickness”  problem  refers  to  the  design  problem  of  infinitely 
thin  blades.  The  blades  are  represented  by  a  distributed  bound  vorticity 
which  is  related  to  the  y-component  of  the  pitch  average  velocity  by  equation 
(2.2-2).  We  seek  the  velocity  field  which  must  satisfy  equation  (2.2-8)  in  the 
cascade  region  and  appropriate  boundary  conditions.  The  blade  camber  lines 
are  solved  iteratively  using  the  blade  boundary  conditions.  The  flow  is 
assumed  to  be  incompressible  and  inviscid.  The  upstream  flow  condition  is 
assumed  to  be  uniform. 


3.1  Flow  regions 

The  flow  field  is  divided  into  three  regions  (figure  1). 

(1)  Region  1 

Region  1  is  defined  in  the  interval  —  oO  <  X  <  0-  In  this  region, 
the  flow  is  irrotational  everywhere.  We  may  write 

=  $(-*)  +  (*/») 

(3.1-1) 

U,T>  r 

where  V  (—  tC  ^  is  the  inlet  flow  condition. 


The  boundary  conditions  in  this  region  are: 


-  far  upstream,  »  O 

-  at  the  leading  edge,  V.  -  V, 

A/  *  aj 

(2)  Region  2 

Region  2  is  defined  in  the  interval  0  <  X  <  1-  In  this 
flow  is  rotational  and  must  satisfy  equation  (2.2*8),  i.e. 

Yi  <-x,  -j)  =  V(*)  + 

£«.  (xlS)  =  V^o)  +  sf*)  V6-(<) 

where  yo<)  is  the  prescribed  pitch  average  velocity. 

The  boundary  conditions  in  this  region  are: 

-  matching  conditions  at  the  leading  and  trailing  edges 

-  no  flow-through  conditions  at  the  blade  surfaces 

(3)  Region  3 

Region  3  is  defined  in  the  interval  1  <  X  <  +  oO  .In 
the  flow  is  irrotational.  We  may  write 

y '>(*•'})  =  y  l+»>)  +  £y  (*,<}) 

-  955,0..,') 

where  ^  ■o')  is  the  outlet  flow  condition. 


region,  the 


(3.1-2) 


this  region, 


(3.1-3) 


The  boundary  condtions  in  this  region  are: 


-  far  downstream,  V/Aj  =  O 

-  at  the  trailing  edge,  \/-  s.  V, 


3.2  Solution  of  the  velocity  field 


In  this  section,  we  present  our  method  of  solving  the  velocity  fields  in 
the  three  regions.  At  the  same  time,  we  try  to  give  reasons  for  our  choice  of 
such  a  method. 

A I  / 

The  solution  for  the  perturbed  velocities  4/^4  are  assumed  to  be  of  the 

form  of  series  of  smooth  functions:  Fourier  series  and  "smoothing”  series 

[Appendix  C].  In  choosing  these  series  solutions,  we  note  that  by  our 

definitions  of  V  ,  equation  (2.2-1),  we  must  satisfy  the  condition  of  zero 

/ \J 

pitch  average  for  A T  .  One  way  of  satisfying  this  condition  is  to  choose 
functions  in  the  series  having  zero  pitch  average.  This  is  one  of  the  property 
of  "smoothing”  functions.  In  addition,  these  "smoothing”  functions  have  other 
properties  which  are  ideal  for  our  applications,  i.e. 


1.  they  are  periodic  in  the  y-direction.  Thus,  they  can  be  used  to 
represent  the  periodicity  of  the  velocity  field  in  the  cascade. 

2.  they  are  derivatives  and  integrals  of  one  another.  This  is  a  useful 
property  for  analysis  purpose. 

IV 

3.  they  are  proportional  to  A  ,  where  n  =  1,2,3,...  .  This 
property  is  desired  in  the  case  of  highly-loaded  blade  where  the 
solidity  is  high  (or  A  is  less  than  1). 

4.  they  have  amplitudes  which  decrease  very  fast  with  increasing 
order  "smoothing"  functions.  Because  of  this  property,  very  few 
terms  in  the  smoothing  series  are  needed  to  represent  a  smooth 
function  while  high  accuracy  can  still  be  achieved. 
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In  the  case  of  infinitely  thin  blades,  we  expand  part  of  the  velocity 

rJ 

potential  0  in  region  2  as  a  series  of  these  "smoothing”  functions.  This  is 
the  ’’smoothing”  technique.  It  turns  out  that  by  using  this  technique,  the 
expressions  for  the  velocity  fields  are  very  simple  to  compute. 

Refer  to  equation  (3.1-2),  we  assume  the  velocity  potential  CD  in  region 
2  of  the  form 

=  &Cx,5)  +  fuCx>9)  (3'lM) 

where  tfis  (X,y)  =  A(y)  1  («)  +  g(<)  T(t)  +  CCX)K(*)  +  ... 

The  coefficients  A (*)  ,  K*)  ,  c.£*)  ,  ...  are  chosen  such  that  the  blade 
boundary  conditions  are  satisfied.  This  form  gives  03  to  have  zero  pitch 
average  and  curl  free. 

The  "homogeneous  part"  of  the  velocity  potential  is  chosen  so  as 

to  satisfy  V  =  0,  to  have  zero  pitch  average,  and  to  match  the 

velocities  at  the  leading  and  trailing  edges  with  those  in  regions  1  and  3. 

The  velocity  field  in  the  cascade  region  can  then  be  written  as 


3.2.1  Smoothing  series 

For  incompressible,  continuity  requires 

V.  V  =  o  <a2-3> 


but 

>. 

X) 

II 

>* 

continuity  of  the  pitch  average  flow),  and 

equation  (3.2-3) 

implies 

V.  £  s 

/v 

O 

(3.2-4) 

or 

= 

-  V.(sV&) 

Now, 

from  equation 

(3.2-1),  we  can  write 

v‘4  = 

-Aiv^r*  so<)  £av«c 

■+■  2  VA .  17*  + 

l__ _ 

+  T(«)[v*A  + 

+  4  75,7*  + 

civ<r] 

+  T  («)  fr's  +  C.  V  « 

• 

• 

• 

+■  2-VC- » ft  -i* 

To  satisfy  (3.2-4),  we  choose 

A)  = 

-  v*.  vs- 

B|Vk  r  = 

-  (Vl6-  +  i7A. 

V*  +  A  lV«cJv 

) 

C  17*)% 

-  (?lA  +  2V6. 

7*  +  5  ft 7«r 

) 

* 

% 

% 

(3.2-51 

and 

the  pattern  for 

all  the  coefficients  is  apparent. 

3.2'.2  Matching  conditions 


TJ . 


As  mentioned  earlier,  the  velocity  potentials  //,  ^  4  are  used  solely 
to  satisfy  the  far  upstream  and  far  downstream  flow  conditions,  and  the 
matching  conditions  of  the  velocities  at  the  leading  and  trailing  edges.  We 
have  chosen  n  ^  A  to  satisfy  Laplace  equation :  therefore 


cf  rr'  -  A*.  I  *  I 

7A<  <•*/!))- L  <•  «• 


we  assume  /^K  A  of  the  forms: 

-UIXl 
<■  « 

1  a.  ( 

CB.«.  +  «,*  )  C  J 

I 

t/  V1  -‘In1) 

\(X,l)=L  Cn  e  e  ^ 

These  assumed  forms  for  7^H  >4  satisfy 


-  Laplace  equation 

-  The  far  upstream  and  far  downstream  flow  conditions 

-  The  periodicity  condition  in  the  y-direction  of  the  flow  field 

The  unknowns  coefficients,  ’  S,  •  B.  and  0 ^  ,  are 

determined  using  the  matching  conditions  of  the  velocities  at  the  leading  and 
trailing  edges: 
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(x  =  o, 

s')  = 

A* 

(*  = 

o,  3 

> 

(X=0, 

a)  = 

% 

t*- 

°l3 

) 

(X*l» 

9)  = 

V 

(x* 

) 

% 

= 

S 

(X- 

'>  3 

) 

Now,  the  velocities  in  the  three  regions  are  the  following. 


(3.2-6) 


For  region  1, 

<?< 

-*> 

II 

V(-co) 

- 

For  region  2, 

il 

K*)  4 

■ 

+  ■St*') 

(V<5-  * 

+■  *(•<) 

(VAr  + 

e  Vie') 

+ 

• 

(vg  + 

CV*  ) 

For  region  3, 

* 

• 

r 

S(— ) 

*  ns 

The  "sawtooth”  function  has  a  finite  jump  across  any  constant 

o(  lines.  At  the  leading  and  trailing  edges,  we  assume  that  the  flow  comes 
and  leaves  smoothly  (zero-incidence  and  Kutta  conditions  respectively).  Thus, 
from  (3.2-7),  we  require  at  X  =  0  and  X  =  1, 
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(  VGr  +  A  V*  )  =  O 


A  (X  *  0,1)  =  O 
Gr'  (X  =  0/  »  O 


(3.2-8) 


From  equation  (3.2-5),  by  choosing  a  loading  distribution  with  the  condition 
Gr  (X  »  Of  l)  =  0,  we  automatically  satisfy  the  condition 

A  (X  a  0,1)  =0. 

Substituting  equation  (3.2-7)  into  equation  (3.2-6)  using  equation  (3.2-8), 
we  can  show  that,  in  region  2 

I =  LL  - T5 - j 

■bso  «ul 

+  (V)  sin  -  (->)*  [A'(0- &«)#<■)]  0D^„  [j 

“  ~  XH  [j-ff+)J -(-»)*  [tfa-CMfk)]  [y 


(3.2-9) 


where 


A,  = 


277 


3.2.3  Blade  boundary  conditions 

The  blade  boundary  conditions  are  used  to  generate  the  blade  camber. 


By  adding  the  blade  boundary  conditions  together,  we  obtain 


Then,  To  represents  the  meam  streamline. 


If  we  define 

the  blade  camber  line  as 

f  ■ 

f.  +  Af 

then,  we  < 

:an  show  that  equation 

(3.2-10)  is  a 

perfect 

derivative  in 

and 

reduces  to 

Af(*) 

=  _!(•)(*!  Af) 

* 

cf) 

—  •  •  • 

-tt 

xso  wsl 

t| 

-  i  [eft) 
L  A*, 

.  W1f(+j]  — 

j  COD  A, 

i  ft)J 

+  ]£*'(*)- w-fki]  +■ 

(d*cw  ... 

A  n 

f(*)J 

(3.2-11) 


3.3  Iteration  process  for  the  blade  camber  line 

We  are  now  in  the  position  to  compute  the  blade  shape  using  equation 
(3.2-11).  We  call  this  equation  the  "camber  generator”.  The  beauty  of  the 


"smoothing”  technique  is  that  the  blade  shape  and  the  velocity  field  can  be 
obtained  by  solving  a  set  of  algebriac  equations,  i.e.  equations  (3.2-5)  and 


(3.2-11): 


A(x)  = 


BOO  * 


Cx)  - 


L 

o*D 


-  G"  +  2. +  A 


('*r> 


-  A*  +  ItW  ■+  B-f" 


(3.3-1) 


-fOO  =  -f.(x)  +  ^f(X) 


(3.3-2) 


Thus,  we  have  (n+1)  equations  and  (n+1)  unknowns.  The  (n+1) 
unknowns  are:  n  coefficients  in  the  smoothing  series  :  Afc)  ,  B(jC)  , 
C  ...  ,  plus  the  blade  camber  line  fCM  . 

We  choose  to  solve  the  blade  camber  line  by  an  iteration  process.  We 
expect  that  the  blade  camber  line  is  not  very  different  from  the  mean 

streamline  -f0  ,  especially  in  the  case  where  the  spacing  <4  is  small.  Thus 
we  can  use  -^0  as  the  initial  guess  for  and  go  through  the  following 

iteration  process: 


1.  Compute  At  ,  B  ,  C  ,  •••  (with  the  guessing  value  of 
Jj.  as  ^  at  the  first  iteration)  using  equations  (3.3-1) 
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2.  Update  -p  by  solving  equation  (3.3-2) 


3.  Go  back  to  step  1  to  update  A 
with  the  new  value  of 
Continue  this  process  until  convergence  in 


JS  ,  C 

computed  in  step  2. 
is  achieved. 

The  whole  flow  field  and  the  pressure  coefficients  at  the  blade  surfaces 
can  be  computed  using  equations  (3.2-7)  and  Bernoulli’s  equations. 


3.4  "Partial  smoothing” 


The  proposed  iteration  process  for  the  camber  line  in  section  3.3  shows 
that  we  can  get  infinite  accuracy  for  -p  by  keeping  an  infinite  number  of 
terms  in  the  smoothing  series.  From  the  engineering  point  of  view,  however, 
we  want  to  take  the  least  number  of  terms  in  the  smoothing  series  for  a  given 
required  accuracy  criteria.  "Partial  smoothing”  refers  to  the  truncation  of  the 
smoothing  series. 


In  the  inverse  poblem,  we  are  most  interested  in  three  quantities:  the 
blade  camber  line  and  the  pressure  coefficients  at  the  blade  upper  and  lower 
surfaces.  Since  -p  +j  V  ,  we  see  that  -p  has  the  same  order  of  accuracy 
as  <^>  .  Similarly,  since  Cp  V  ,  we  see  that  Cp  also  has  the  same 
order  of  accuracy  as  •  We  now  construct  a  table  showing  the  truncation 
errors  in  -p  and  C-  as  a  function  of  the  number  of  terms  used  in  the 
smoothing  senes  <p  . 


Number  of  terms  kept 
in  the  smoothing  series 


Truncation 

errors 


*(«) 

*(o) 


3.5  Numerical  method 

As  an  example  of  the  design  problem,  we  write  a  program  to  compute 
the  blade  camber  line  by  keeping  two  terms  in  the  smoothing  series.  In  this 
case,  the  smoothing  series  is  of  the  form: 

@  (*/•?}  *  *(*)  T(*)  +  &C<)3~(«') 

Let’s  investigate  the  order  of  magnitude  of  the  truncation  error  in 
-j-~  when  the  above  "partial  smoothing”  form  is  used.  From  equation 
(3.2-11),  we  see  that  the  truncation  error  in  is  of  the  order  of  K  (o)  , 

or  (Appendix  C|: 


error  in 


•foo 


7ZO 


Therefore,  for  the  above  example,  the  truncation  error  in  ~f~  is  about  three 

_ i _  _ l _  .1 _  X  V"  D..  1 _ : _ _ 1-.  , _  . _ • 


orders  of  magnitude  less  than 


By  keeping  only  two  terms  in 


the  smoothing  series,  we  can  still  maintain  very  high  accuracy  (accurate  to 
.OOltf  A*  )•  Moreover,  the  numerical  calculation  is  very  simple.  Given  the 
loading  distribution  (x)  as  an  analytical  form,  the  blade  camber  can  be 
calculated  through  a  set  of  algebriac  equation.  No  numerical  technique  is 
needed  to  compute  derivatives  or  integrals. 
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The  iteration  process  for 
parameters: 


i 


involves  calculating  the  following 


1.  Guess 


f 


and  its  derivatives  as: 


-foo  =  +  L  (g'-ai') 

* 

jU  =  -f.'oo 
j‘(*)  -  f.\*) 


for  simplicity,  we  set 
iteration  process 


H  // 

(x)  =  £  (*.)  through  out  the 


2.  Compute 

A(X)  = 


*'(*)  = 


a’*' 

o+rr 


B(<) 

3.  Update 


-  6"  +  z-f  V  +  A-f 

o+r) 


j'(x)  =  -fo  (*)  +  i  (5-"-  A-f'.  Aj')  -d(FsOM) 
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f  00  =  -f«0O  +  ji-  ( 6- A -f')  -  FS0*l 


where 


F5U 


'  •£.  -XhIX-*M 


M  =  ZZ  ■S"T» - •  (-')V)c«Xn[-f(^-fO)] 

■fcsO  nr/  ^*1 

i 


4.  Go  back  to  step  2  and  repeat  the  calculations  until  convergence  in 
-j.  is  achieved. 

The  convergence  criteria  used  in  the  calculation  is,  at  each  location 

»!♦  I  **» 

i  <*0  -  f(*)  <  ERROR 

i.e.  if  the  difference  between  the  present  value  of  and  the  previous  value 

of  at  all  locations  X  i  is  less  than  ERROR  (the  convergence  criteria), 

then  convergence  in  is  achieved. 

3.4  Numerical  results 


Numerical  results  shows  that  comvergence  in  can  be  achieved  fast, 

depending  on  the  spacing  to  chord  ratio  A  .  In  this  section,  we  present  a 
numerical  example.  The  results  are  compared  with  the  "exact”  solution  (using 
the  Biot-Savart  law)  (4j. 


We  have  taken  the  case  of  an  inlet  guide  vane  with  the  following  inputs: 

-  spacing  to  chord  ratio  A  =  0.75 

-  inlet  angle  =0° 

-  outlet  angle  —  45° 

-  loading  distribution  4P  X(l-X) 

Calculations  are  made  at  21  points  and  the  convergence  criteria  ERROR 
is  JO 

17  iterations  are  required  for  -f-  to  converge.  The  computer  program 
output  is  shown  in  table  1.  Figure  2  show  the  plots  of  the  blade  camber  line 
and  the  mean  streamline.  Note  that  there  is  reverse  curvature  of  the  blade 
near  the  trailing  edge.  This  observation  is  also  found  in  the  results  performed 
by  Tan,  Wang,  Hawthorne  and  McCune  [4].  The  computational  time  for  this 
example  is  around  7  seconds  CPU  time  on  the  Digital  VAX-11  computer. 

Table  2  shows  a  comparason  of  the  blade  camber  obtained  from  the 
"smoothing”  technique  and  the  "exact”  method.  The  results  compare  well. 

This  preliminary  study  of  the  "smoothing”  technique  shows  that  it  is 
indeed  a  very  powerful  engineering  tool  for  the  design  problem  in  terms  of 
computational  time  and  accuracy  (if  desired). 
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Chapter  4  :  The  "Finite-Thickness”  Problem 


In  this  chapter,  we  will  present  a  method  to  solve  the  inverse  problem 
for  blades  having  finite  thickness.  The  method  will  be  similar  to  the  one 
presented  in  the  ”  Zero-Thickness”  problem,  except  that  we  now  have  to 
differentiate  the  average  flow  quantities.  In  the  case  of  blades  having  finite 
thickness,  we  can  define  two  average  quantities: 


1.  The  pitch  average  velocity  which  is  defined  as: 


The  pitch  average  velocity  has  no  physical  representation  of  the 
flow  since  there  is  no  flow  in  the  blade.  However,  if  we  are  to 
model  the  presence  of  the  blades  by  distributing  singularities  on 
the  blade  camber  lines,  then  V  does  exist  in  this  sense. 


2.  The  gap  average  velocity  which  is  defined  as: 

•f+4-T 

—  /  1  1 
vT  (x)  = 

A#  • 


r f+i“T 

/  Jvi 

Jci-r  J 


4  -  2T(tc)  +  t 

Obviously,  the  gap  average  velocity  does  represent  the  average  of 
the  actual  flow. 


In  the  design  problem  of  blades  having  finite  thickness,  as  mentioned  in 
the  introduction,  we  choose  to  model  the  presence  of  the  blades  by  distributing 

vortices,  sources  and  sinks  on  the  blade  camber  lines.  In  our  method,  the  two 

—  / 

main  given  quantities  are:  the  loading  distribution  VT^  ,  and  the  blockage 
distribution  T" 

Given  these  two  quantities,  the  gap  average  velocity  is  known.  It  is 


defined  by: 


where  (x)  is  found  using  continuity,  i.e. 

4 

4  -  AT(jf) 

By  distributing  vortices  on  the  blade  camber  lines,  we  have  shown  that 
the  strength  of  the  vorticity  is  related  to  the  y-component  of  the  pitch 
average  velocity  through  equation  (2.2-2).  Therefore,  in  the  "Finite-Thickness” 
problem,  we  no  longer  know  the  strength  of  the  vorticity.  Equation  (2.2-8)  is 
still  valid  in  describing  the  flow  field  in  region  2,  however,  we  now  choose  a 
different  way  to  describe  it. 


^Tx  M 


4.1  Flow  regions 

Using  the  same  general  approach  as  in  the  ” Zero-Thickness”  problem,  we 
again  divide  the  flow  field  into  three  regions.  The  velocity  potentials 
k  ^  are  use(*  to  satisfy  ^ar  upstream  and  far  downstream  flow 
conditions,  and  the  matching  conditions  of  the  velocities  at  the  leading  and 
trailing  edges.  We  again  make  use  of  the  "smoothing”  functions  to  satisfy  the 
blade  boundary  conditions. 

The  three  flow  regions  are: 

(1)  Region  1 

Region  1  is  defined  in  the  interval  —  <  X  <  0.  In  this  region, 


the  flow  is  irrotational  everywhere.  We  may  write 
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y,  =  yc--) + 

«,»)  =  ^(X,,)  (4-l-D 

where  y(~>  is  the  inlet  flow  condition.  The  boundary  conditions  are: 

-  far  upstream  -  O 

-  at  the  leading  edge  V  4  -  V 

(2)  Region  2 

Region  2  is  defined  in  the  interval  0  <  yC  <  1.  We  choose  to  analyze 
the  flow  in  a  region  between  the  blade  camber  lines.  The  flow  in  this  region 
is  divergence  free  and  curl  free.  We  may  write 

7<z5  =  o  <41-2> 

The  boundary  conditions  in  this  region  are: 

-  matching  conditions  at  the  leading  and  trailing  edges 

-  no  flow-through  condition  at  the  blade  surfaces 


(3)  Region  3 

Region  3  is  defined  in  the  interval  1  <  X  <  +  *0  In  this  region, 
the  flow  is  irrotational.  We  may  write 

y>ov?>  =  y<>~)  +  £>(*>•}) 

jfjOOs)  = 

where  is  the  oulet  flow  condition.  The  boundary  conditions  in  this 


region  are: 
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far  downstream 


O 


-  at  the  trailing  edge  V  s  V. 


4.2  Solution  of  the  velocity  field 


Consider  the  flow  between  the  two  camber  lines  (excluding  these  camber 

p  +  _ 

lines)  located  at  u  =  -j-  and  Lj  ==  A  A  .  In  this  region,  we  can 


write  as  a  velocity  potential  satisfying  Laplace  equation: 

yu *>*)  * 

VV  :  O  I4-2'1 1 

Let’s  assume  to  be  of  the  form 

=  «5(^  GyC*)  +  ?>£*) 

+  rfrjCXfy)  +  9)  (4-2‘2* 

where  the  smoothing  series  ^rs  is  assumed  to  be  of  the  form 

0TS(*,y)  *  &rsCx)  + 

4  At  (x)  X  (<)  +  4  -  *  *  (4.2-3) 

and  thus,  0  can  be  written  as 

$(*/*£)  ~  [&T  *•“  "$>  +  ^T5  * 

at  r(*)  +  sT  cr$0  4  --  <4-2*4) 

Combining  equation  (4.2-1)  ,  (4.2-2)  along  with  the  choice  vX  -  0. 
we  obtain 


V  s  -  +  2  V  <?y .  ^7*.  -  GtVJS  (4.2-5) 


Based  on  the  experience  from  the  "Zero-Thickness”  problem,  there  are  a 
few  motivations  in  choosing  ^  of  the  above  form: 


in  the  smoothing  series  A.  is  the  additional  amount  of 
vorticity  needed  to  represent  the  total  amount  of  vorticity 
distributed  on  the  blade  camber  lines. 


and  represent  some  average  in  the  potential  velocity 

which  we  have  the  freedom  to  choose  at  our  convenience. 


is  chosen  to  satisfy  the  matching  conditions  at  the  leading 
and  trailing  edges. 


-  Ar  ,  ST  ,  CT  .  —  are  chosen  to  satisfy  the  blade  boundary 
conditions. 


We  are  now  in  the  position  to  solve  for  the  velocity  field  in  region  2. 
Our  tasks  are  to: 


1.  relate  some  of  the  unknowns.  in  the  assumed  form  of 
given  gap  average  velocity  VT 

2.  satisfy  the  Poisson  equation  (4.2-5) 

3.  satisfy  the  blade  boundary  conditions 

y,  ■  =  o 

yt  .  vc*+t)  =  o 

4.2.1  Relation  between  the  gap  average  velocity  and  the 

unknowns 

In  this  subsection,  we  relate  some  of  the  unknowns  in  the  velocity 
potential  d)  defined  in  equation  (4.2-4)  to  the  gap  average  velocity. 


<&  to  the 


(4.2-6) 

(4.2-7) 


By  definition  of  the  gap  average  velocity, 


i 


I  Xt  ^  *  U-J. r)VT 

7*t  - 

=  2  SCT )  VT 

Evaluating  the  above  integral  using  equation  (4.2-4),  we  obtain 

C  &  +j '«r  )  **  -  C-frr)  ia  +  (V«ST4*V^)T  = 

We  now  choose  the  coefficients  in  ^  such  that 

=  1T 

Therefore,  we  require 

r)  =  VVx 


(  ^  )  =  O 


Equation  (4.2-8)  requires 


&'  *  VTX  -  -f^T 


Equation  (4.2-9),  a  vector  quantity,  can  be  satisfied  by  choosing 


Define 


■»X 


•) 


(4.2-13) 


We  will  show  later  that 


relates  to  the  source/sink  strength. 


4.2.2  Satisfaction  of  the  Poisson  equation 

We  now  choose  the  coefficients  in  the  smoothing  series  so  that  the 
Poisson  equation  (4.2-5)  can  be  satisfied. 

From  the  assumed  form  of  the  smoothing  series  ,  equation  (4.2-3), 

we  can  write 

+  s(*)  (vV  +  at  vh  4-  2  VA7.  ft  *  gT  /Pit I*  ) 

*  !(*)  (VlAr  4-  6,  V*  +  J.  VgT .  Br  +  CT  /P*!'’  ) 

*■  TW  (7 lBr  +  Cr  V*  *■  /  flCr.Wr  +J>r  IfcT) 

■t  .  .  . 


Substituting  the  above  equation  into  equation  (4.2-5)  using  equations  (4.2-10) 
through  (4.2-12),  we  obtain,  in  scalar  form 


Ar  IMV  -  (  <x+/)  -  f(*T*t  ) 

-  S(*)  [  (G-t  *'*')-  At  f  -  lArf  *  &r  IV*  !* ] 
4  I(»)  [  Ar  -  2<f'+  Crl?*r] 

*  T(«)  [  6'V  -  CTf-  J>r  /*/’] 


The  above  equation  is  of  the  form 

*.(*)  -  b‘i'))a-4Cx)  *  ^tCx)  4- 

One  way  to  satisfy  the  above  equation  is  to  choose 

«.<*)  *  o 

fti  (*)  ~  ° 

c,t  (*)  =  O 


Using  the  choices  suggested  above,  we  obtain 

>  . «-»  it.  ^  ^  «»  /  \ 


At  |P*  |V  -  f  (G-r  +  S  )  +  (Vrx  +/) 


= 

CT  /(fc/1*  = 
J>r  /P*/*  = 


-  Ar  +  5rf  + 

-  5'V  +-  cr  f  "  * 


/  */ 


//  _  /  / 
4-  <£  0-r  j- 


(4.2-14) 

By  choosing  AT  ,  UT  ,  CT  ,  ...  of  the  above  forms,  it  turns  out  that 
one  of  the  blade  boundary  conditions  is  satisfied.  This  fact  will  be  discussed 
in  section  (4.2.4). 


4.2.3  Matching  conditions 


v. 


.4s  in  the  infinitely  thin  blade  case,  the  velocity  potential  7 far)  A  are 
used  to  satisfy  the  matching  conditions  of  the  velocities  at  the  leading  and 
trailing  edges,  and  the  far  upstream  and  far  downstream  flow  conditions. 
Again,  are  chosen  to  satisfy  Laplace  equation.  A  similar  analysis 

for  as  in  the  "Zero-Thickness”  problem  gives: 


U>  V  V'  -  4 *  I*-*  I 

L  (*/•))  =  Z-  Z-  73 

trO  n*\  'S 


tsO  n*| 


.2-15 


With  the  analysis  done  so  far,  we  can  now  write  down  the  velocity  in 


region  2: 


y*.  (*,*)  =  yr(x)  +  [^X)ix 

+  5  (x )  +  $(x ))  -•-  Ar(x  )  Vot  J 

+  T(pt)  [Vat(x)  +  bT(x)  Vot  J 
+  T  (•*]  £v&r(*)  (-7 '{*)  V*  _J  +■  ‘ 

(4.2-16) 


Again,  we  require 

At  *  O 

/ 


at 


X  =  0,1 

£-T'+S'  =  O 

for  the  flow  to  come  and  leave  smoothly  leading  and  trailing  edges  (the  zero- 
incidence  and  the  Kutta  conditions  respectively).  Using  the  definition  of 
A-y-  in  equation  (4.2-14),  the  above  conditions  become: 


/ 


—  / 


VT*  *fi 


/ 


-  o 
*  o 


at  X  =  0,1 


(4.2-17) 


4.2.4  Blade  boundary  conditions 

So  far,  we  have  chosen  all  the  coefficients  in  the  expression  for 
We  must  now  satisfy  the  two  blade  boundary  conditions,  i.e.  equations  (4.2-6) 
and  (4.2-7).  It  turns  out  that  one  boundary  condition  is  automatically  satisfied 
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through  the  choices  of  Ar  >  &T  »  Cr  ,  ...  in  equation  (4.2-14)  while  the 
other  blade  boundary  condition  is  used  to  generate  the  blade  camber  line. 

By  adding  and  subtracting  equations  (4.2-6)  and  (4.2-7),  the  blade 
boundary  conditions  become: 

<  >T  .  V«.  -  A t(70).VT  *  o 

At  (V0) .  >T  .  pr  =  o 

Evaluating  at  the  upper  and  lower  sufaces  of  the  blade  and 

substituting  them  into  the  above  equations,  we  obtain: 

-C&r+l)  +  (?&+  V?n  )  ■  >r  ■  V* 

+  X  Cr)  [vAt-V«  *  Br/M1] 

+  K  (t)  £  VCr  .  *  -&r  /  ®  I  J  ■+■  ... 

=  At  .  VT 

+  «5CT)  ^  •  P^"  1 

L  J  (4.2-18) 

■+-  T(t)  »  P*T  -r  CTP*.  P7”  J  •  •  • 


and 


Ar(nj.  V* 

+  SCr )  [VC 6-T  +  thVx  +  ArIV«r] 

+  *tr)  fveT,&<  +-  cT  /r«r]  +■ 

=  -  (g-t^)v«.vt  +(v0+v?n).vr 

+  <  >T  .  VT 

+  1(T)  [<7Ar  .  VT  +  0r  Wf .  Pr  J  (4.2-19) 

+  k(t)  [vcr,vr  +  3>rv«  -VT ]  +  . .. 

By  substituting  the  choices  of  4<y  i  5y-  ,  Cr  •••  *n  equations  (4.2-14) 
into  equation  (4.2-19)  and  along  with  the  identity 

A  “  At(7/aJ  .  ttc  -  <  7  it>r.VT 

we  can  show  that  the  blade  boundary  condition  (4.2-19)  is  automatically 
satisfied. 

Finally,  the  other  blade  boundary  condition,  equation  (4.2-18),  is  satisfied 
by  choosing  the  blade  camber  line  ~f  appropriately.  We  will  call  equation 
(4.2-18)  the  "camber  generator”. 

Substituting  the  choices  of  the  coefficients  Af  ,  dT  ,  Cr  ,  ...  into 
the  blade  boundary  conditions  (4.2-18),  we  obtain 


scr)VT.V«  =  'V* 

*  scr)  |AT(7^).Vr  -  <  y£;> .  ft 

4  Hi  4^i  "s/} '  M'  I 

♦[*«*  *  ^  <«J{  *  -  Crf '] 

*  *  “  •  (4.2-20) 

This  is  a  differential  equation  relating  -p  and  the  coefficients  Ay- , 
5r  ,  C  x  >  •••  •  As  in  the  "Zero-Thickness”  problem,  we  can  reduce  the 
above  equation  into  an  algebriac  equation  for  -p  (in  the  iterative  sense). 

Again,  define 


Then  yo  represents  the  mean  streamline. 

Also,  define 

•foo  =  -foOO  ■*  Af(*) 

Then,  referring  to  equation  (4.2-20),  it  can  be  shown  that 

^Ct)Vt  .  =  -  r  i  (Af> 


[c(r)d  +  yjL  T(r) 

=  -  [?(r)+sMT(r)][(<sS+S')-4rf' J 

J 

[LCT)di  *  y£<(-)]{6'r-crrj 

»  i.  [j.(r)+Sfr>K(T')jr  Br  -  CTf] 

<*K 

* 


and  finally 


r 

•SO-)'  At  (v£ ) .  yr  -  <  >T  .  7*  ■ 
c  - 

-  XT)  C^jsTV  •  vr  =  dL  (raw) 


where 


„  . .  v1  v'  <-~x'  |,ut|  f  *>» 

F50H  5  2-  Z-  75 - (  — 7 -  “**  SCt)oosA„T'  ), 

■fcaO  nsl  ***  \  An  y 

[  (-!)%(*)  -  1  [a;  (t)  .  Cr(i)f'(t)] 

-  ^  3T(*)  +  . . .  ]  cos  A.  [fw-  f(*)1 

Ai|  J 

[[Ar(*)-6r(t)|'(t)J  1-  ^  cr(t) 

A  * 

Combining  these  equations  together,  the  '"camber  generator”  becomes 


*5lh  An  ^ 


-f-  SCt)  cos  AmT" 


■fM  =  -  £  JV(t)  +  seO  Kt')]  (e-r+t'-Arf') 

-  j  [UT)  +  5(r)  K(r)l(g'  -Ct{')  . 


~  (  FSOv)^) 
A 


(4.2-21 


4.3  Iteration  process  for  the  blade  camber  line 


Oi 


We  choose  a  very  similar  iteration  scheme  used  in  the  ”  Zero-Thickness 

problem.  The  iteration  process  consists  of  calculating  the  coefficients  A~j~ 

Cy  / 

3t  ,  CT  ,  ..  using  equations  (4.2-14).  Then,  update  b  and  usin 

equations  (4.2-11)  and  (4.2-13).  Finally,  update  -j-  using  equation  (4.2-21). 


In  quantitative  forms, 


1.  Compute 


A'(G-'W)  +  (v' */) 


aTf  -  ^Brf 


^  +f) 


-tyx)  = 


.  -A-  1  Crf  +  ICrf 


0+f,v) 


2.  Update 


(fi 


(4.3-2) 


-f(X)  =  'fc(X')  +  Af(*) 

In  order  to  start  the  above  iteration  process,  we  need  to  know  not  only 
an  initial  guess  for  ,  but  also  for  ^  and  %  A  discussion  of  the 

"physical”  reprsentations  of  ^  and  T  is  given  in  Appendix  D. 

4.4  Numerical  method 

A  computer  program  is  written  to  solve  the  inverse  problem  using 
"partial  smoothing”  (refer  to  section  3.4).  As  mentioned  in  section  4.3,  the 
iteration  process  for  -j-  consists  of,  in  "partial  smopthing”  forms: 


1.  Guess 


/■  f'-  t 


.  .  .  .*  .  * 


//  -f 

p  =  VTX  (  sin  2irX  -  l) 

s'  =  o 

f  *  -fo(*1  +  J  £t(t)  -f5(T)X(T)J  (G-'+i'-Arf') 


2.  Compute 


^t(*)  = 


(<rr  -f  S  )  *  (  vy,  +/Q 

U+f1') 


5  Oc)  ,  7&.+y)  ±  ±dl  1  lh 
r( ’  (I*. fl) 

3.  Update 


/  .  f 


f(V)  =  -  L.  £t(t)  +  +^-  \j') 

-  i: 


4.  Check  for  convergence.  If  not,  go  back  to  step  2  to  recompute 
4r  ,  with  the  updated  values  of  ,  "%/  and 

and  continue  this  process  until  convergence  in  -L.  is 
achieved.  ’ 

As  in  the  "Zero-Thickness”  problem,  the  convergence  criteria  for 
~p  used  during  the  iteration  process  is,  at  every  location  X-’L  ,  we 


require 


n+i  vt 

f  (*-‘)  -  -f  (*0  < 


ERROR 


for  convergence  to  succeed. 

After  convergence  has  been  achieved,  we  proceed  to  calculate  the  pressure 


coefficients  on  the  blade  surfaces  defined  by: 


where  V  is  defined  in  equation  (4.2-16),  and  *onst-b  is  defined  as 


o*»s  cl 


k  = 


I  t 

I  +■  -jp  4- 


Finally,  in  order  to  accelerate  convergence,  we  pre-calculate  all  leading 
and  trailing  edges  variables  needed  in  the  iteration  process.  The  flow  chart 
for  the  computer  program  is  shown  in  table  4. 

The  above  iteration  process  requires  computing  derivatives.  Two  methods 
of  computing  derivatives  were  investigated:  the  Spectral  method  (Chebyshev 
collocation)  [7],  [8],  and  the  finite  difference  method  (central  difference).  The 
finite  difference  method  was  chosen  over  the  Spectral  method  method  because 
it  is  numerically  more  stable. 

Numerical  problems  were  encountered  in  the  iteration  process  for  -j-  . 
Due  to  "partial  smoothing”,  this  method  is  not  able  to  resolve  the  singular 
point  in  the  source/sink  distribution  at  the  leading  edge.  It  was  found  that 
when  approximately  11  points  (depending  on  the  value  of  the  spacing  A  ) 
are  used  in  the  calculation,  convergence  in  -£■  is  achieved  in  about  10 
iterations.  When  more  than  11  points  are  used  in  the  calculation,  the 
iteration  process  fails  to  converge.  We  propose  to  use  a  filter  to  resolve  this 
problem.  Studies  of  this  numerical  problem  and  the  filtering  method  are 
discussed  in  Appendix  E. 


4.5  Design  choices 


In  our  method,  there  are  two  main  input  parameters  available  to  the 
designer:  The  blockage  distribution  /  )  and  the  loading  distribution 


VT  (X)  •  As  an  example,  we  use  analytical  forms  for  both  and 

V^(x)  33  inputs  to  our  numerical  code. 


4.5.1  Blockage  distribution 

In  our  numerical  example,  the  blockage  distribution  T(x)  is  chosen  to 
be  of  the  form 

T(x)  ^ 


Oj  k 

X  (/-x) 


We  define  the  maximum  blockage  parameter  BLOCK  as 
BLOCK  = 


^  r 7 

We  restrict  ourselves  to  the  case  where  /  =  0  at  the  trailing  edge. 

If  T  ( 4  )  ^0,  then  a  stagnation  point  must  exist  at  the  trailing  edge. 
We  do  not  think  that  this  is  a  good  model  of  the  real  flow.  In  the  actual 
low  speed  flow  situation,  the  potential  flow  outside  the  boundary  layer  is 
pushed  away  from  the  trailing  edge  by  the  presence  of  the  wake.  Therefore, 

we  think  that  the  condition  T/(4)  =  0  is  a  better  model  for  the  real 

flow  giving  more  realistic  pressure  coefficients  at  the  trailing  edge.  This 

condition  require  cu  >1.  We  note  that  this  is  not  the  restriction  of  our 

method,  b  can  be  any  real  number. 


4.5.2  Loading  distribution 


In  our  numerical  example,  the  loading  distribution  is  taken  to  be  of  the 

form 

s  -  ap  -  x 

where  c,  and  cL  can  be  any  real  numbers. 

With  these  design  choices,  we  will  be  able  to  study  some  effects  of  the 
blockage  and  loading  distributions  on  the  pressure  coefficients  at  the  blade 
surfaces.  A  more  practical  way  would  be  to  read  in  the  thickness  and  loading 
distributions  at  discrete  points  and  use  a  numerical  method  to  compute  their 
derivatives  and  integrals. 


4.0  Numerical  results 


In  this  section,  we  will  first  discuss  the  limitations  of  our  current 
numerical  code.  Then,  we  will  attempt  to  close  the  loop  for  our  method  using 
a  direct  method.  Finally,  some  results  are  presented. 


Limitations  of  the  current  numerical  code 

Our  current  numerical  code  (using  the  iteration  scheme  of  section  4.4)  is 
limited  to  the  following  cases: 


1.  when  the  axial  chord  is  divided  into  10  intervals  or  less  (<  11 
points),  -p  converges  without  using  the  filter.  Otherwise,  the 
filter  is  needed  for  convergence  in  ~p  to  succeed. 


2.  when  the  spacing  to  chord  ratio  A  is  of  the  order  of  1  and 
greater,  -f  converges  slow  and  fails  to  converge  when  the 
convergence  criteria  ERROR  is  less  than  /O”*’. 

These  problems  can  be  resolved  if  many  more  terms  in  the  smoothing 
series  (equation(4.2-3))  are  kept  so  that  the  singular  point  at  the  leading  edge 
can  be  resolved.  However,  we  decide  not  to  do  so  because  the  idea  of  the 
’’smoothing”  technique  is  to  be  able  to  achieve  high  accuracy  using  very  few 
terms  in  the  smoothing  series. 

Closing  the  loop 

Our  design  method  is  supposed  to  find  the  blade  shape  which  is  supposed 
to  do  two  specified  jobs:  a  certain  amount  of  circulation,  and  a  certain 
loading  distribution.  Given  these  two  parameters  and  the  blockage 
distribution,  our  numerical  code  computes  the  corresponding  blade  shape  and 
pressure  coefficients  at  the  blade  surfaces. 

Does  this  blade  profile  actually  do  the  specified  jobs?  In  an  attempt  to 
answer  this  question,  we  use  a  direct  method  to  compute  the  circulation  and 
the  pressure  coefficients  of  the  blade  shape  obtained  from  our  indirect  method. 
If  there  are  agreements  in  these  results  between  the  two  methods,  then  we 
succeed  to  close  the  loop  for  our  design  method. 

We  choose  the  direct  method  developed  by  McFarland  [10].  It  makes  use 
of  the  panel  method.  The  numerical  code  has  many  options,  some  of  which 
are  what  we  need  to  perform  the  comparason. 

In  general,  results  of  the  comparason  show  that: 

-  for  the  individual  pressure  coefficients,  there  are  good  agreements 

between  the  two  methods  (  ^  h%)  away  from  the  leading  and 


trailing  edges  (.1<  X  <.9).  Near  the  edges,  the  two  results  differ 
substantially  (see  example  discussed  below). 

-  for  the  loading  (or  pressure  difference  across  the  blade)  distribution, 
the  two  results  agree  within  5%. 

-  for  the  circulation  (or  comparing  the  oulet  angle),  the  two  results 
agree  well  within  5%. 

An  example  of  closing  the  loop  for  our  design  problem  is  now  presented 
We  have  taken  the  case  where: 

O 

-  inlet  angle  CXj  =  0 

-  outlet  angle  =45° 

-  BLOCK  =  .1 

-  spacing  /4  =  .  5 

-  t  ^  x  0  _x)v 

-  Ap  ~  X°'S(/-x) 

Table  4  and  5  show  the  numerical  results  of  our  indirect  method  using 
11  points  and  41  points  respectively.  Note  that  the  two  results  agree  well 
within  the  error  of  the  numerical  scheme  used  to  compute  derivatives.  Figures 
3  and  4  show  the  corresponding  blade  shapes  and  pressure  coefficients. 

The  pressure  coefficients  for  the  above  blade  shape  (obtained  from  our 
indirect  method)  are  calculated  using  the  direct  method.  Two  options  in  the 
direct  method  program  are  used: 

1.  Option  1  -  the  inlet  flow  condition  is  specified,  and  the  program 
uses  the  Kutta  condition  to  determine  the  outlet  flow  condition. 

2.  Option  2  -  both  the  inlet  and  outlet  flow  conditions  are  specified 
(or  the  circulation  is  specified). 


Table  6  shows  the  results  obtained  using  option  1.  The  results  show 

that: 

-  the  exit  flow  angle  is  45.42  ,  compared  to  the  specified  outlet  angle 

of  45°. 

-  the  loading  distribution  is  of  the  specified  shape  (figure  5). 

Table  7  shows  the  results  obtained  using  option  2.  The  results  are  very 
similar  to  those  obtained  using  option  1. 

Figure  6  shows  a  comparason  of  the  pressure  coefficients  obtained  from 
the  direct  and  indirect  methods.  It  shows  that  the  two  results  agree  well 
away  from  the  leading  and  trailing  edges.  Near  these  edges,  the  pressure 
coefficients  obtained  from  the  direct  method  show  oscillations.  At  the  trailing 
edge,  because  our  blade  shape  is  thin,  the  direct  method  fails  to  fit  a  curve 
through  the  control  points  giving  negative  thickness.  Consequently,  the  Kutta 
condition  at  the  trailing  edge  is  not  satisfied,  and  the  results  near  the  edges 
are  not  to  be  trusted. 

We  conclude  that  the  blade  shape  obtained  from  our  design  method 
succeed  to  do  the  two  specified  jobs.  We  were  unable  to  close  the  loop 
completely,  but  the  individual  pressure  coefficients  obtained  from  the  two 
methods  do  have  the  same  general  shape. 

Comparason  of  the  "Zero-Thickness"  and  "Finite-Thickness” 
results  in  the  zero  blockage  limit 

It  is  found  that,  in  the  limit  of  T (x')  going  to  zero,  the  blade  camber 
lines  obtained  from  the  two  theories  are  not  the  same.  Although  the  two 
theories  themselves  are  the  same  in  this  limit,  the  results  are  not  the  same 
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because  the  "Zero-Thickness”  problem  does  not  predict  the  presence  of  a 
stagnation  point  at  the  leading  edge  while  the  "Finite-Thickness”  problem  does. 

Table  8  shows  an  example  of  how  the  camber  lines  obtained  from  the 
two  theories  compare.  Plots  of  these  two  blade  camber  lines  are  shown  in 
figure  7. 

Effects  of  the  spacing  to  chord  ratio  on  the  blade  camber 

Figure  8  shows  the  effects  of  the  spacing  to  chord  ratio  on  the 

blade  camber  line.  As  A  increases,  a  smaller  number  of  blades  are 
available  to  do  the  same  job  resulting  in  highly  cambered  blades. 

Effects  of  the  maximum  blockage  on  the  pressure  coefficients 

Figure  9,10,11  show  the  effects  of  increasing  the  parameter  BLOCK  (the 
maximum  blockage)  on  the  pressure  coefficients.  We  have  taken  the  case 
where: 

-  inlet  angle  =  O 

_  a 

-  outlet  angle  =  ^  5 

-  spacing  to  chord  A  —  0,5 

.  T  ~  XO-X')'' 

-  Ap  -  X(i-X) 

As  the  parameter  BLOCK  increases,  the  flow  near  the  maximum 
thickness  location  (in  this  example  X  =  .33333)  accelerates  due  to  the 
Venturi  effect.  This  effect  can  result  in  highly  unfavorable  pressure  gradients 
at  both  the  upper  and  lower  surfaces  if  the  loading  distribution  is  not  chosen 


properly.  Figure  11  shows  such  a  situation.  A  better  loading 
distributionshould  be  sought  for  this  example  to  give  more  acceptable  pressure 
coefficients. 


Rounded  leading  edge  example 

Figure  12  shows  an  example  of  an  inlet  guide  vane  having  a  rounded 
leading  edge.  Compare  with  figure  10  (this  is  a  fair  comparason  since  all 
inputs  in  the  two  cases  are  the  same  except  for  the  blockage  distribution),  we 
see  that  the  effect  of  the  rounded  leading  edge  is  to  further  overexpand  the 
fluid  near  the  leading  edge.  However,  we  do  not  think  the  current  numerical 
code  can  resolve  the  rounded  leading  edge  case  accurately  because  of  the 
limitations  on  the  size  of  the  computational  intervals. 


4.7  Design  procedure 


The  numerical  code  outlined  in  section  4.5  can  be  used  to  explore  the 
effects  of  blade  loading  and  blockage  distributions  on  the  behavior  of  the 
pressure  coefficients  at  the  blade  surfaces.  Given  some  design  requirements, 
the  designer  can  use  the  above  program  to  find  the  "best”  blade  shape  using  a 
trial  and  error  method. 

In  this  section,  we  present  an  example  of  a  design  procedure  using  our 
numerical  code.  Suppose  that  we  wish  to  design  inlet  guide  vanes  which  can 
do  the  following  jobs: 


o  • 

-  the  flow  is  to  be  turned  from  ctfj  =  0  to  =  45  ,  with  a 

spacing  to  chord  ratio  A  =  .5. 


-  the  blockage  distribution  is  to  have  an  analytical  shape  of  the  form 
X  (  I  —  X  )*”  with  the  maximum  blockage  parameter 
BLOCK  =  0.1. 

We  seek  a  loading  distribution  which  gives  "good”  pressure  coefficients  at 
the  blade  surfaces  (in  terms  of  minimizing  flow  separation).  We  start  the  trial 
and  error  process  by  choosing  the  loading  distributions  which  is 

1.  highly  loaded  near  the  leading  edge  (maximum  at  =  .33333). 

Figure  13  shows  the  corresponding  blade  shape  and  pressure 

coefficients. 

2.  highly  loaded  at  midchord.  Figure  14  shows  the  corresponding 
blade  shape  and  pressure  coefficients. 

3.  highly  loaded  near  the  trailing  edge  (maximum  at  X*.  ==  .66666). 

Figure  15  shows  the  corresponding  blade  shape  and  pressure 

coefficients. 

Comparing  the  results,  we  conclude  that,  for  the  above  example,  case  1 
gives  the  "best”  pressure  coefficients.  In  general,  results  show  that  by  loading 
high  near  the  leading  edge,  the  corresponding  pressure  coefficients  are  "good”. 
Figure  16  shows  an  example  of  a  compressor  blade  after  going  through  the 
above  ♦rial  and  error  process.  Figure  17  shows  an  example  of  an  impulse 
blade  going  after  going  through  the  same  process. 


Chapter  5  :  Conclusion 

A  two-dimensional  design  method  for  highly-loaded  blades  was  presented 
in  this  thesis.  Singularities  are  distributed  on  the  blade  camber  lines  to  model 
the  presence  of  the  blades.  The  non-linear  theory  is  based  in  part  on  the 
Clebsh  formulation.  A  "smoothing”  technique  was  used  to  solve  for  the  blade 
boundary  conditions.  Numerical  examples  was  presented  using  a  "partial 
smoothing”  form  of  the  iteration  scheme  for  the  blade  camber  lines. 

It  was  found  that  when  the  blades  are  assumed  to  be  infinitely  thin,  the 
blade  camber  lines  can  be  solved  through  an  iteration  process  of  a  set  of 
algebriac  equations.  The  iteration  process  converges  very  fast  seconds 

CPU  time  on  the  Digital  VAX-11  computer)  for  the  typical  solidity  range 
found  in  turbomachines.  The  results  compare  very  well  with  those  obtained 
from  an  "exact”  method. 

When  the  blades  are  assumed  to  have  finite  thickness,  the  "partial 
smoothing”  form  of  the  iteration  scheme  for  the  blade  camber  lines  fails  to 
resolve  the  singular  point  at  the  blade  leading  edge  accurately.  In  order  to 
get  high  accuracy,  an  infinite  number  of  terms  in  the  smoothing  series  would 
have  to  be  kept,  making  the  "smoothing”  technique  less  attractive  compare  to 
other  techniques.  A  practical  numerical  code  based  on  a  "partial  smoothing” 
form  of  the  iteration  scheme  for  the  blade  camber  lines  was  presented  giving 
very  fast  convergence  solutions  (<^10  seconds  CPU  time  on  the  Digital 
VAX-11  computer)  with  satisfactory  accuracy. 
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INLET  ANCLE  2  0  000 
OUTLET  ANCLE  2  45  000 
NUNS Eft  OF  POINTS  2  21 

PARABOLIC  LOADING  INPUT  PROPORTIONAL  TO  x(l-x) 


CONVERGENCE  HISTORY  OF  FIX) 


ITER  *  1  -  £RRMAX  *  0  01456  AT  X  2  0  90000 

ITER  *  2  -  ERRNAX  2  0  00336  AT  X  2  0  60000 

ITER  *  3  -  ERRNAX  *  0  00127  AT  X  2  0  45000 

ITER  •  4  -  ERRNAX  *  0  00039  AT  X  *  0  55000 

ITER  »  5 - ERRNAX  2  0  00027  AT  X  2  0  45000 

ITER  »  6  -  ERRNAX  2  0  00013  AT  X  a  0  40000 

ITER  »  7  -  ERRNAX  =  0  00011  AT  X  2  0  45000 

ITER  *  8  -  ERRNAX  *  0  00007  AT  X  =  0  SOOOO 

ITER  »  9  -  ERRNAX  2  0  00006  AT  X  *  0  45000 

ITER  »10  -  ERRNAX  *  0  00004  AT  X  2  0  SOOOO 

ITER  *11  -  ERRNAX  *  0  00003  AT  X  2  0  45000 

ITER  *12  -  ERRNAX  2  0.00003  AT  X  2  0  SOOOO 

ITER  *13  -  ERRNAX  2  0  00002  AT  X  2  0  45000 

ITER  *14  -  ERRNAX  2  0  00002  AT  X  2  0  SOOOO 

ITER  *15  -  ERRNAX  2  0  00001  AT  X  2  0  45000 

ITER  *16  -  ERRNAX  2  0  00001  AT  X  2  0  SOOOO 

ITER  *'7  -  ERRNAX  2  0  00001  AT  X  2  O.SSOOO 


X 

FN(X) 

FLOW  ANCLE 

FIX) 

BLADE  ANCLE 

0  00000 

0  00000 

0  00000 

0  00000 

-9  10520 

0  05000 

0  00012 

0  41538 

•0  00809 

-9  27549 

0  10000 

0  00095 

1  60386 

-0  01633 

-8  61382 

0  1S000 

0  00312 

3  47844 

-0  02324 

-6  47547 

0  20000 

0  00720 

S  93741 

-0  02768 

-2  95814 

0  25000 

0  01367 

8  88065 

•0  02841 

2.22304 

0  30000 

0  02295 

12  18862 

•0  02380 

8  79703 

0  35000 

0  03537 

15  73517 

-0.01293 

IS  43656 

0  40000 

0  05120 

19  39206 

0  00382 

21  19087 

0  45000 

0  07062 

23  03761 

0  02583 

26  04913 

0  SOOOO 

0  09375 

26  56504 

0  05269 

30  19280 

0  S5000 

0  12062 

29  88813 

0  08403 

33  76094 

0  SOOOO 

0  15120 

32  94323 

0  11954 

36  83603 

0  65000 

0  18537 

35  68779 

0  15893 

39  48203 

0  70000 

0  22295 

38  09645 

0  20192 

41.75705 

0  75000 

0  28367 

40  15600 

0  24821 

43  68834 

0  80000 

0  30720 

41  36034 

0  29744 

45  27878 

0  85000 

0  35312 

43  20571 

0  34918 

46  50883 

0  90000 

0  40095 

44  18653 

0  40285 

47  31881 

0  05000 

0  45012 

44  79155 

0  45763 

47  54620 

1  00000 

0  SOOOO 

45  00000 

0  5121 6 

47  42213 

Table  1  Numerical  example  of  ’Zero-Thickness’  problem 


64 


INPUT 

SPACING  S  *  0  75 

INLET  ANCLE  «  0  000 

OUTLET  ANGLE  »  45  OOO 

PAAA80LIC  LOADING  INPUT  PROPCHTIOfAL  TO  x(l-x) 


X 

EXACT 

F 

SMOOTH 

F 

EXACT 

BLADE  ANGLE 

SMOOTH 
BLADE  ANGLE 

0  00000 

0  00000 

0  00000 

-6  36362 

-9  10520 

0  05000 

-0  00733 

-0  00809 

-8.77228 

-9  27549 

0  10000 

-0  01498 

-0  01633 

-7  83866 

-8  61382 

0  15000 

-0  02102 

-0.02324 

-S  31557 

-6  47547 

0  20000 

-0  02428 

-0  02768 

-1  64080 

-2  95814 

0  25000 

-0  02390 

-0  02841 

2  90811 

2  22304 

0  30000 

-0  01923 

-0  02380 

8  06026 

8  79703 

0  35000 

-0  00979 

-0  01293 

13  51490 

15  43656 

0  40000 

0  00475 

0  00382 

18  96298 

21  19087 

0  45000 

0  02450 

0  02583 

24  13298 

26  04913 

0  50000 

0  04947 

0  05269 

28  83133 

30  19280 

0  55000 

0  07947 

0  08403 

32  95549 

33  76094 

0  SOOOO 

0  11423 

0  11954 

36  48344 

36  83603 

0  £5000 

0  15338 

0  15893 

39  44676 

39  48203 

0  70000 

0  19650 

0  20192 

41  90393 

41  75705 

0  75000 

0  24316 

0  24821 

43  91742 

43  68834 

o  aoooo 

0  29288 

0  29744 

45  53462 

45.27878 

0  85000 

0  34517 

0  34918 

46  77274 

46  50883 

0  90000 

0  39940 

0  40285 

47  59450 

47  31881 

0  95000 

0  45473 

0  4S763 

47  84251 

47  54620 

1  00000 

0  50954 

0  51216 

46  68512 

47  42213 

Title  2 


Comparison  of  'smoothing'  and  ’exact*  results 


Inputs 


Table  3 


Floe  chart  for  the  'Finite-Thickness'  proble* 


INPUT  PARAMETERS 


UAX  BLOCKAGE  =  3  10003 

SPACING  =  0  S00C0 

INLET  ANCLE  =  0  00000 

OUTLET  ANCLE  =  45  00000 

NUMBER  OF  POINTS  IJX  =  11 

MAX  NUMBER  OF  ITERATIONS  ALLOWED  =  10 

MAX  ERROR  IN  F(X)  ALLOWED  EIIR11AX  =  0  000001 

FILTERING  OPTION  *  0 

BLOCKAGE  AND  LOADING  PARAMETERS 
A  =  1  00 

B  -  2  00 

C  =  0  SO 

D  *  1  00 


ITERATION 

f 

i  — 

- ERRUAX 

=000376 

AT 

X 

at 

0 

10000 

ITERATION 

• 

2— 

- ERRUAX 

=0  00048 

AT 

X 

m 

0 

10000 

ITERATION 

• 

3— 

—  ERRUAX 

=0  00013 

AT 

X 

3 

0 

10000 

ITERATION 

• 

4  — 

—ERRUAX 

<0  00002 

AT 

X 

a 

0 

10000 

ITERATION 

1 

5-  — 

—  ERRllAX 

=0.00000 

AT 

X 

s 

0 

10000 

ITERATION 

< 

5--- 

- ERRUAX 

=0  00000 

AT  X 

» 

0 

10000 

X  LOADCX)  T(X)  Fll(X)  F(X)  Cp»  Cp- 


0 

00000 

0 

00000 

0 

00000 

0 

00000 

0 

00000 

1 

00000 

1 

00000 

0 

10000 

1 

09742 

0 

01367 

0 

00259 

-0 

01371 

0 

21609 

-0 

83669 

0 

20000 

l 

37954 

0 

02160 

0 

0133S 

-0 

00982 

0 

32123 

-1 

00828 

0 

30000 

1 

47839 

0 

02481 

0 

03600 

0 

01488 

0 

35421 

-0 

91697 

0 

4000C 

1 

46323 

0 

02430 

0 

07164 

0 

05381 

0 

30674 

-0" 

E6721 

0 

£0000 

1 

36328 

0 

02109 

0 

12073 

0 

10615 

0 

21504 

-0 

85185 

0 

60000 

1 

19472 

0 

01620 

0 

18316 

0 

17181 

0 

09294 

-0 

83185 

0 

70000 

0 

96733 

0 

01063 

0 

25817 

0 

24986 

-0 

05235 

-0 

80053 

0 

80000 

0 

68977 

0 

00540 

0 

34412 

0 

33873 

-0 

21089 

-0 

74507 

0 

30000 

0 

36581 

0 

00152 

0. 

43835 

0 

43595 

-0 

37577 

-0 

66037 

1 

00000 

0 

00000 

0 

00000 

0 

53713 

0 

53751 

-0 

55867 

-0 

55867 

Table  4 


Kuoencal  example  of  'Finite-Thickness’  problem 
(without  filter) 


INPUT  PARAMETERS 


MAX.  BLOCKAGE  •  0.10000 

SPACING  ■  0.50000 

INLET  ANGLE  ■  0.00000 

OUTLET  ANGLE  •  45.00000 

NUMBER  OF  POINTS  I JK  *  41 

NA„'  NUMBER  OF  ITERATIONS  ALLOWED  *  10 

MAX.  ERROR  IN  FIX)  ALLOWED  ERRMAX  •  0.000100 

FILTERING  OPTION  •  1 

BLOCKAGE  AND  LOADING  PARAMETERS 
A  •  1.00 

B  •  2.00 

C  •  0.S0 

D  •  1.00 


ITERATION  # 

1 

- ERRMAX  *0.01190 

AT  X  • 

0.02500 

ITERATION  » 

2 

- ERRMAX  *0.00676 

AT  X  • 

0.02500 

TERATION  i 

3 

- ERRMAX  *0.00388 

AT  X  * 

0.02500 

ITERATION  * 

4 

- ERRMAX  *0.00223 

AT  X  • 

0.02500 

ITERATION  # 

5 - ERRMAX  =0.00127 

AT  X  * 

O.O2S0O 

ITERATION  * 

6 

- ERRMAX  *0.00074 

AT  X  • 

0.02500 

ITERATION  * 

7 

- ERRMAX  *0.00041 

AT  X  • 

0.02500 

ITERATION  » 

8 

- ERRMAX  *0.00024 

AT  X  • 

0.02500 

ITERATION  * 

9- 

- ERRMAX  *0.00013 

AT  X  ■ 

0 . 02500 

ITERATION  * 

10 

- ERRMAX  *0.00009 

AT  X  • 

0.02500 

X 

LOAD  IX) 

T  (X) 

FMlX) 

FIX) 

Cp* 

Cp- 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

1.00000 

0.02500 

0.58000 

0.00401 

0.00009 

-0.00398 

0.50002 

0.24987 

0.05000 

0.79921 

0.00761 

0.00048 

-0.00772 

0.00085 

-0.61389 

0 . 07500 

0.95307 

0.01083 

0.00133 

-0.01066 

0.03782 

-0.66715 

0. 10000 

1 .07077 

0.01367 

0.00273 

-0.01245 

0 . 08002 

-0.71882 

0. 12500 

l . 16390 

0.0161S 

0.00474 

-0.01297 

0. 13205 

-0.75553 

0. 15000 

1.23856 

0.01829 

0.00741 

-0.01217 

0. 17935 

-0.79703 

0. 17500 

1.29845 

0.02010 

0.01076 

-0.01010 

0.22816 

-0.83309 

0.20000 

1 . 34604 

0.02160 

0.01484 

-0.00680 

0.27418 

-0.86304 

0.22SOO 

1 . 38307 

0 . 02280 

0.01966 

-0.00234 

0.31633 

-0.S8633 

0.254100 

1.41086 

0.02373 

0.02523 

0.00322 

0.34S97 

-0.90261 

0.27500 

1 . 43040 

0.02439 

0.03159 

0.041982 

0.354ia 

-0.91 166 

0.30000 

1.44248 

0.02481 

0.03873 

0.01742 

0.34936 

-0.91 355 

0.32500 

1 . 44776 

0.02499 

0.04667 

0.02599 

0.32412 

-0.90880 

0.35000 

1 . 44677 

0 . 02495 

0.05543 

0.03550 

0.309O1 

-0.89827 

O.37SO0 

1 . 43995 

0.02472 

0.06501 

0.04592 

0.30709 

-0 . 88303 

0.40000 

1 . 42769 

O.O2430 

0.07541 

0.05724 

0.29402 

-0.86644 

0.42500 

1.41031 

0.02371 

0  08664 

0 . 06943 

0.27274 

-0.85046 

0.45000 

1 .38810 

0.02297 

0 . 09870 

0.08249 

0.24809 

-0.84213 

0.47500 

1.36131 

0.02209 

0. I l 159 

0.09639 

0.22259 

-0.3388 S 

0.50000 

1 .33017 

0.02109 

0. 12531 

0.  1  11  12 

0. 19656 

-0.83721 

0.52500 

1 .29487 

0.01999 

0. 13986 

0. 12667 

0.  16941 

-0 . 83899 

0.55000 

1.255S8 

0.01879 

0. 15523 

0. 14303 

0. 14059 

-0.S3399 

0.57500 

1.21248 

0.01753 

'  0. 17141 

0.  16016 

0.  10993 

-0.82397 

0.60000 

1 . 16570 

0.01620 

0. 18838 

0.  17804 

0.07753 

-0.82351 

0.62500 

1 . 11538 

0.01483 

0.20613 

0. 19667 

0 . 04360 

-0.81760 

0.65000 

1.06164 

0.01344 

0.22465 

0.21600 

0.00838 

-0.81086 

0.67500 

1 .00-159 

O.U12U3 

0  24390 

0.23602 

-0.02793 

-0 . 80297 

0.70000 

0.94453 

0. n 1063 

0  26387 

0.25670 

-0.06516 

-0.79369 

0 . 72500 

0 .  SKi  r  if, 

II. l>4  >9  ">5 

0  28452 

0.27800 

-0. 10315 

-0.78285 

0.7SOOO 

0. SI  456 

(1.1*179) 

0  30582 

0.29990 

-0. 14174 

-0.77030 

0.775«  m 

t>.  74522 

n  00662 

0.32773 

0.32237 

-0. 18080 

-0 . 75593 

O.nuomo 

0.67302 

0.00540 

0 . 35021 

0.34537 

-0.23020 

-0.73968 

0.82S00 

0 . 59802 

0. 004  26 

0,37320 

0 . 36888 

-0.25988 

-0.721  S3 

O.SSOOO 

0 . 52030 

0.00323 

0,39667 

0.39287 

-0 . 29980 

-0.70153 

0.87500 

0 . 4399 1 

0.00231 

0.42055 

0.41729 

-0.34003 

-0 . 67980 

0.90000 

0.35692 

0.00152 

0.44478 

0.44213 

-0. 38069 

-0 . 65653 

0.92SOO 

0.27138 

0 . 00088 

0.46931 

0.46732 

-0.42202 

-0.63202 

0.9SOOO 

O. 18335 

0.00040 

0.49406 

0.49282 

-0.46446 

-0.60669 

0.97500 

0.09237 

O.UOOlO 

0.51896 

0.S185S 

-0.50899 

-0.58137 

1 .0004)0 

0.004X10 

0.04X100 

0.S4394 

0.54434 

-0.55891 

-0.SS891 

Table  5  Numerical  example  of  ’Finite-Thickness'  problem 
(vith  filter) 


CONTROL 

CONTROL  POINT 

COORDINATE 

SOURCE 

NORMAL 

TANGENTIAL 

PRESSURE 

POINT 

X 

Y 

DLNSITY 

VELOCITY 

VELOCITY 

COEFFICIENT 

1 

0.987487 

0.531406 

-0.327305 

0.000000 

-1.003151 

-0.006311 

2 

0.962497 

0.505435 

-0.329590 

0.000000 

- 1 . 579879 

-1  .496010 

3 

0.937506 

0.479427 

-0.328757 

0.Ol)»**)0 

- 1 . 1 736  n 

-0.377364 

4 

0.912511 

0.4S351S 

-0.326162 

O.OOOMO0 

-1.3179U3 

-0.736.867 

5 

0.S875I5 

0 . 427783 

-0.322227 

0.00**100 

- 1 . 3(14845 

-0.702619 

6 

0.862519 

0 . 402295 

-0.317136 

0.004)000 

-1 .312181 

-0.721820 

7 

0.837522 

0.377107 

-0.310993 

o .  coonoe 

-1.316798 

-0.733958 

8 

0.812525 

0.352270 

-0.30384)9 

0.0*  X*M  )0 

-1 .321*121 

-0.745096 

9 

0.787528 

0.327831 

-0.295573 

0 . 000**10 

-1.324733 

-0.754919 

10 

0.762531 

0.303837 

-0.286240 

0.000000 

-l .328028 

-0.763659 

1 1 

0.737534 

0.280335 

-0.27576O 

0  .  0O*MK)0 

-l .330985 

-0.771520 

12 

0.712537 

0.257369 

-0.2644)55 

0.000000 

-I.3337o2 

-0.778762 

13 

0 .687539 

0.234986 

-0.2SIO53 

0.000000 

-1.336256 

-0. 785580 

1 4 

0  662542 

0.213232 

-0.236684 

0.000000 

-1.338725 

-0.792184 

IS 

0.637544 

0.  192150 

-0.220876 

0.000000 

-1.341169 

-0.798733 

16 

0.612S4? 

0. 171783 

-0. 203552 

0.004**00 

- 1 . 343636 

-0.805356 

17 

0.587549 

0.  152174 

-0. 184650 

0.00**100 

-1.346141 

-0.81 2096 

18 

0.562550 

0. 133362 

-0. 1641 10 

0.000***) 

-1 .348704 

-0.8I9OO2 

19 

0. 537552 

0.  1  15385 

-0. 141866 

0.000000 

-l .351329 

-0.826091 

20 

0.512553 

0.098279 

-0. 117855 

0.00**100 

-l  .3540*)2 

-0.833322 

21 

0.487553 

0.082078 

-0.092026 

0 . 000*100 

- 1 . 356702 

-0 . 840639 

22 

0.462554 

0 .  *1668  15 

-0.064315 

0.  OO**)O0 

-1.359402 

-0.84797S 

23 

0.437553 

0.OS2S20 

-0.034663 

0.0000*10 

-l .362056 

-0 . 855 1 98 

24 

0.412553 

0.039225 

-0 . 003003 

0.o*H*mo 

-  1 . 36460 1 

-0.862137 

25 

0  387552 

0.026961 

0.030739 

0.000000 

-  1 . 366939 

-0868 523 

26 

0.362550 

0.015759 

0.066643 

0.000(100 

- 1 . 368936 

-0.873986 

27 

0.337548 

0.005653 

0. 104813 

0.000**10 

-l.37o39t 

-0.877972 

28 

0  312545 

-0 . 003320 

0.  14536 1 

0.0**1000 

-1 .371003 

-0.879651 

29 

0.287541 

-0.01 1 1 18 

0. 188412 

0.000000 

- 1 . 37*1336 

-0.877822 

30 

0.262536 

-0.01 7693 

0.234096 

0.000000 

- 1 . 367744 

-0.870723 

31 

0.237530 

-0 . 0229*90 

0.282508 

0. OOOO00 

- 1 . 362293 

-0.8S5842 

32 

02 12523 

-0.026944 

0.333655 

0.00*1000 

-1.352654 

-0.829674 

33 

0. 187514 

-0.029482 

0.38735 l 

0 . 00*1000 

-1.337004 

-0.787580 

34 

0. 162502 

-0.030524 

0.443059 

0.000000 

-1.31 2949 

-0.723835 

35 

0. 137490 

-0.02999 2 

0 . 499690 

0.000000 

-1.277559 

-0.632158 

36 

0. 1 12476 

-0.027823 

0. 555407 

0  00*3*100 

-1 .22 ‘>688 

-0.507219 

37 

0 . 087465 

-0.023995 

0.607521 

0 . 0*10*100 

-1 . 159363 

-0.344(23 

30 

0 . 06246  l 

-0  OI9S70 

0.6S2563 

0.0***  *00 

-1.073582 

-0. 152578 

39 

0.O37470 

-0.011764 

0.685564 

0  01***10 

-0.908008 

0. 17S357 

40 

0.01 4985 

-0.004*66 

0 . 704)727 

0.01**100 

-0.993459 

0.013040 

4| 

0  002499 

-0.000828 

0.709200 

0. 001**10 

-0.201021 

0 . 95959 | 

42 

O.0O25O0 

0.  i*NM)l  1 

-0 . 4-48950 

0.******) 

O.SK7S20 

0.654821 

43 

0  OIS0H0 

0. 0*10041 

-0.4SI883 

0.0****10 

0.925071 

0.  144244 

44 

0 . O37500 

-0.000047 

-0.457399 

0.OO***H> 

0.765250 

0.414393 

4S 

0  062501 

-0.00*8)39 

-0. 442719 

0.OOO**» 

0.771431 

0.44)4895 

46 

0.087505 

0.000588 

-0.414927 

0  01***10 

0.7S*696 

0.42X923 

47 

0. 1 1 25  10 

0.iX»2OH| 

-0.381 196 

0.Ol*MMtO 

0.745731 

0 . 443X85 

40 

0.  137514 

6.004531 

-0.345486 

0  ******) 

0  741692 

0 . 449S93 

49 

0.  162518 

0 .  *107946 

-0.310129 

0.OOO1*)0 

0.  742463 

0  448749 

50 

0. 187521 

0.01 2296 

-0  276276 

0.  U*NM*10 

0  746.854 

0.442210 

51 

0.212523 

0.017535 

-0.244315 

0.0****)O 

0.753821 

0.431754 

52 

0.237525 

0.023614 

-0.214207 

0.Oo***)O 

0.762583 

0.418467 

53 

0.262526 

0.030492 

-0.  18S719 

0.000000 

0.772603 

0.403084 

S4 

0.287527 

0.038138 

-0. 158562 

0.000000 

0.783S36 

0 . 3X607 1 

55 

0.312529 

0.046526 

-0. 132463 

0.01**100 

0 . 70S  1 73 

9.367700 

56 

0.337530 

0 . 05564 1 

-0. 107195 

0.00**100 

0.807401 

0.348104 

57 

0.362531 

0.065472 

-0.082593 

0.0*10000 

0.820164 

0.327331 

58 

0.387533 

0 . 0760 1 S 

-0. 058547 

0. 0**1000 

0.833439 

0 . 305379 

59 

0.412534 

0.087267 

-0.034989 

0.0**1000 

0.847222 

0.282214 

60 

0.437536 

0.099238 

-0.01 190! 

0  .  *)u*  *  too 

0.661516 

0.257791 

61 

0.462537 

0.111 899 

0.010708 

0 . 0****)0 

0.876309 

9.233082 

62 

0.487538 

0.  125281 

0.032818 

0.000000 

0.891593 

0  205062 

63 

0.512539 

0. 139374 

0.0S4374 

0 .  O***H)0 

0.907342 

0 .  1 7673 1 

64 

O.537S40 

0.154177 

0.075340 

0. 0*10*100 

0.923510 

0. 147129 

65 

0.562540 

0.  169688 

0.095655 

0  O*M*IO0 

0.94«*)SI 

0. 1 16304 

66 

0.587541 

0.  1859*)3 

0.  1  15273 

0 . 00*  *  *  10 

0.956890 

0.084361 

67 

0  612541 

0.202814 

0 . 1 34 l S7 

0.00*1000 

0.973947 

0.051427 

68 

0.637541 

0.220415 

0. 152273 

0. 00**100 

0.991135 

0.0 17650 

69 

0.662540 

0 . 238695 

0. 169601 

0.OO**)O0 

l . 0083S8 

-0.016786 

70 

0.687S40 

0.257642 

0.  186(36 

0.000000 

1  .025518 

-0.051686 

71 

0.712539 

0.277243 

0.201884 

0. 0****)0 

1.042525 

-O . 086858 

72 

0.737539 

0.2974*6 

0.216868 

0 . 0****)0 

1 . 059300 

-0.122116 

73 

0.762538 

0.3 I83S6 

0.231 109 

0.0****)0 

1 .075786 

-0.  1573  IS 

74 

0.787537 

0.339838 

0.244641 

0.  0****l0 

1.091 943 

-0. 192339 

75 

0.812536 

0.361918 

0.257491 

O.O****)0 

1 . 107753 

-0.2271 18 

76 

0.937535 

0 . 38458 1 

0.269679 

0. 0»***H) 

l .  123289 

-0.261 779 

77 

0.862534 

0.40781 1 

0.281 197 

0. l)**K)O0 

l . 1 3XO90 

-0.295250 

78 

0.887533 

0.431S89 

0292010 

0  00**100 

1  .  156313 

-0.3370S9 

79 

0.912531 

0. 45589 1 

0.302015 

0.00**) 00 

l .  150741 

-0.324205 

80 

0.937527 

0.484)694 

0.31 1029 

0  ****100 

t . 227682 

-0.597203 

81 

0.962521 

0.505916 

0.31863S 

0.000**10 

1  .  189794 

-0.415610 

82 

0.987512 

0.531485 

0  323538 

0.0UUUO0 

1 .242S29 

-0.543878 

UN  I  POM  EURO*  (EPS)  •  -« .  747S7S4B-03 


CT.  •  -«.90»W2S6 

CD  •  0  2866O34E-02 

CIRCULATION  ■  -0.4525191 

POTENT  I AL  FLOW  VELOC1TT  DIAGRAM 

UPSTREAM  VEI.OCITY  ■  0.T9I75  AT  0.00004  DEGREES 

ONSET  VEIOCITT  •  1.00000  AT  28.905-46  DEGREES 

DOWNSTREAM  VELOCITY  •  1.27056  AT  45.J2360  DEGREES 

NOTE:  ALL  VELOCITY  QUANITIES  ARE  SCALED  BY  THE  ONSET  VELOCITY  ••• 


Table  6 


Results  froa  direct  Method 
(inlet  condition  specified) 


69 


CONTROL 

CONTROL  POINT 

COORDINATE 

SOtR'E 

SOR8M 

TANGENTIAL 

PRESSURE 

POINT 

N 

V 

DENSITY 

VI  1.0*  1  TY 

VELOCITY 

cum  IC1ENT 

1 

0  '3874x7 

0. 511406 

-0  332912 

0  IMNMM40 

-1 .01 1283 

-0 . 022692 

2 

0.962497 

0.505435 

-0.335192 

*4  (NHNH>0 

-l  .52*»S96 

-1 .31221 l 

3 

0.937506 

0.479427 

-0.334161 

0 .0*4*4*  »K> 

-  1 .  1X0087 

-0.392604 

4 

0.31251 1 

0.453515 

-*>.331771 

0  *4lR4*W40 

-1 .34*6018 

-0.705682 

5 

0.NH7SIS 

0.4277S3 

-0.327845 

O.iNHHNH* 

- 1 . 297752 

-0  6S4 161 

6 

0. 862519 

0.4O229S 

-0. 3-2764 

0  IMMlIH) 

-1.3'  >638  1 

-0.706631 

7 

U. 837522 

0.377107 

-0316634 

0  INNNM40 

-l  312422 

-0.722451 

8 

0.812525 

0.352270 

v  -4)  34*9463 

O.IMMHMH* 

-1 .317907 

-0. 7368x0 

9 

0. 787528 

0.327831 

-0.3**1243 

1*  (NMNNK) 

- 1 . 322764 

-0.749704 

10 

0. 762531 

0.3*4.1837 

-0.291927 

*>  IMMMMM) 

-1 . 327<*x5 

-0.761 153 

1  1 

0.717534 

0.2*41335 

-4*  28  1  465 

0 . 1 N  N  M  i|  *0 

-l  .33**945 

-0.771416 

12 

0.712537 

0  257169 

-0  269780 

0 

-1.334444 

-0. 7SO740 

13 

0  687539 

0. 2349^6 

-4*  256799 

|>  (NMHMH) 

-1.337655 

-0  7S9322 

14 

0.662542 

0.213232 

-0.2 42451 

O.lMNNIiM) 

-1  . 34i  *666 

-0.797387 

*5 

0.637544 

0.  192154* 

-♦) .  226666 

o  inmhmm) 

-  t  .  34354.8 

-0.M4.S123 

16 

0.612547 

0.  |7Pn3 

-0.2*4' >16  5 

I).(NMMNN) 

-1 .346364 

-0.812696 

17 

0.587549 

0.152174 

-0.  194*4x5 

0.O4HN1440 

-1.349143 

-0.82*4187 

IN 

0.562550 

0.  133162 

-4*.  ;  69968 

o .  ihhmnH) 

-1 .351919 

-0.82-685 

19 

0.537552 

0.  1  15385 

-O.  14  7  744 

0  (M  M  N  H  H) 

-1 .3547 io 

-0.8.15239 

20 

0. St  2553 

0  149X3?9 

-** .  123  753 

o  U*KK»K> 

-  1 . 3575 l 1 

-4.842x37 

21 

0.487553 

0 . 41x20  78 

-0  097940 

0  INN  M  NH) 

-1 .36**310 

-*)  850443 

0.462554 

0  4if.6x IS 

-♦*  i|7ii242 

1)  iknMN) 

-i .363088 

-0  858*  *08 

23 

0  437553 

0.052520 

-44  IMitfJNl 

0  .  (JiNMNN) 

-  1 . 36580 1 

-4.865413 

24 

0.442553 

0.4*39225 

-0.»H  18944 

1)  IKMIMI 

-1 .36S391 

-0.872495 

25 

O  387552 

0.026  361 

0.O248INI 

o  .  |N  N  MM  M* 

-l  .37**762 

-0.87X990 

26 

0  362554) 

0  015759 

0.  *>6*17  |  4 

0  .  IM  MM  M  Ml 

-1 .372782 

-0.884530 

27 

0.337548 

0 . 008653 

*»  «*989»i2 

O.OIMHNN) 

-1.374249 

-0.8SX560 

2* 

0.312545 

-0.0*41320 

O  139479 

0  (MM  MMH* 

-1.374x63 

-0 . 89* » 247 

29 

0  287541 

-0  0  1  l  l  1  8 

O  ls<2S74 

0  INMNMM) 

-1  374184 

-0.888J82 

30 

0.262536 

-0.01 7693 

**  228315 

0  1MMMH40 

-1 .371565 

-0.881191 

31 

0.237530 

-0.022990 

0.2'68«>3 

0  IXMMM40 

- 1 . 366069 

-0 . 866 1 44 

32 

0.212523 

-0.026J44 

0  32m >48 

0.  IHMMM40 

- 1 . 356364) 

-0.839713 

33 

0.  187514 

-0 . 029482 

0.381x65 

0.  INMNMM* 

- 1 . 340608 

-0.797229 

34 

0.  1 6254 >2 

-0 . O3o524 

0.437723 

0  1  M  N  M  M  *0 

-1 .3164**9 

-0.732933 

35 

0  1374  JO 

-0.029992 

0.494534 

0  IMMMMM) 

- 1 . 2S0822 

-0.640505 

36 

0. 1 12476 

-O  027823 

0  S5*4455 

O.IMMMMM) 

-  1 .23**687 

-0.51 4590 

37 

0.4187465 

-0.023995 

0  602789 

0  inmnmH) 

-1  .  162'MJO 

-0.35**243 

38 

0.062461 

-0.01 8570 

0  648*448 

0.04MMH10 

-1.075764 

-0. 157269 

39 

O.H37470 

-0.011  764 

0  6X1224 

0  (*l  M  N  M*0 

-0.9*  >9  196 

0. 173363 

40 

0.014985 

-0 . 0*14866 

0  696473 

0  0*  M  M  M*0 

-0.995337 

0.009303 

41 

O  4103499 

-0  44008  28 

0  704997 

0 . 0*  •*  M  M  lO 

-0. 192771 

0962839 

42 

0.4N4254 40 

0.0**4  401  | 

-0.443632 

0.IMHHM4I* 

0  59  54402 

0.645973 

43 

0  0  1 

0.0**4  4*141 

-0.446573 

0  IMMMMM) 

0.929794 

0. 135483 

44 

0. 03754  40 

-0  0*40047 

-0.4521446 

0.  O0*MMM* 

0. 768674 

0.409140 

45 

0.446254  4  1 

-O.OU*H)39 

-0.437382 

0  .  IMMNNIO 

0.774482 

0 . 40*11 78 

46 

0.  ‘7875*45 

0  4)00588 

-0.4*19512 

0  O*NNN40 

0.758529 

0.424634 

47 

0. 112510 

0.41**24  48 1 

-0.375695 

O.INHM10O 

0.748417 

0.439872 

48 

0. I 375 1 4 

0  ****4531 

-0 . 3399**4 

0  .  U*  M  M4*)0 

0.744294 

0.446042 

49 

0. 162518 

0 . 007946 

-O.  3**4474 

0  .  *****  M  4l*0 

0.744995 

0.444982 

50 

0. 187521 

0.012296 

-0.270560 

O.OI40UOO 

0.749350 

0.438474 

51 

0.212523 
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Table  7 


Results  from  direct  method 
(circulation  specified) 


IMPUT 


BLOCK  *  0  OOOOl 

SPACING  S  =  0  75 

INLET  ANCLE  =  0  900 

OUTLET  ANGLE  =  46  000 

PARABOLIC  LOASINC  INPUT  PROPORTIONAL  TO  »(!-«) 


X 

FIX) 

finite 

thicks*** 

remit 

FIX) 

zero 

thicks*** 

remit 

0  00000 

0  00000 

0  00000 

0  06000 

-0  00849 

-0  00809 

0  10000 

-0  01803 

-0  01633 

0  15000 

-0.02680 

-0  02324 

o  aoooo 

-0  03334 

-0  02768 

0  25000 

-0  03620 

-0  02841 

0  30000 

-0  03347 

-0  02380 

0  3S000 

-0  03376 

-0  01293 

0  40000 

-0  00761 

0  00382 

0  45000 

0  01410 

0  02583 

0  soooo 

0  04083 

0  05269 

0  55000 

0  07215 

0  08403 

0  60000 

0  10770 

0  11951 

0  65000 

0  14715 

0  15893 

0  70000 

0  19020 

0  20192 

0  75000 

0  23654 

0  24821 

0  80000 

0.28581 

0  29744 

0  85000 

0  33757 

0  34918 

0  90000 

0  39123 

0  40285 

0  96000 

0  44600 

0  45763 

1  00000 

0  50053 

0  51216 

Table  8  :  Comparason  of  'Zero-Thickness*  and  'Finite-Thickness' 
results  in  the  zero  blockage  liiit 
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Figure  2  Blade  camber  and  lean  streamline 
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Figure  3  Blade  shapes  obtained  from 
’Smoothing’  technique 
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Figure  8 


Effects  of  spacing  to  chord  ratio 
on  the  blade  casbers 
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Figure  9  Effects  of  aaxiauB  blockage  on  Cp's 
(BLOCK  =  0  001) 
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Figure  10  Effects  of  aaxiaia  blockage  on  Cp's 
(BLOCK  =  01) 
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Figure  17  :  Impulse  blade 
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Appendix  A  :  The  Relationship  between  the  Swirl  Schedule  and 
the  Pressure  Difference  across  the  Blade 

Consider  the  cascade  geometry  shown  below 


Under  the  assumption  of  incompressible,  inviscid  and  uniform  inlet  flow 
condition,  Bernoulli’s  equation  is  valid  everywhere.  We  can  write: 

p+  +  •  r  *  i 

Therefore 

AP  *  ?*-?'  x  ( \T+  i/+)(  V'.  V+)  (A-i) 

We  can  do  the  following  approximation: 

Cv-  +  ^  z  vT 


To  estimate  (  \/  ri  we  consider  the  circulation  along  the  closed 

path  0  -  (2)  -  (5)  -  @  shown  in  the  figure.  From  Kelvin’s  theorem,  we  can  write 


Now, 


z  •  n  +  ft  +  ps  +  p,  =  o 

i*4.T 

Ja  s  (-1'  *^T)  VT„  (X<4X) 

V  J 

r?  S  -  V". 

4-  V  A/ 

PJ  *  -  ^(V/D)  <=>3  z-(4-*T)V '  (*) 

^4T  * 

I7,  =  V'+.  dl*  «  V*  AX 


where 


nr  i 


AX-Jlf  (^'iT')*  *  AX 


(A-2) 


Substitute  the  above  relations  into  (A-2),  we  obtain 

CV-.  V*  )  AX  x  (4-iT)[vT^X*M).  VT;|t<)] 

and  thus,  in  the  limit  of  AX  _ p  O 

cv-v-)  *  vT;. 

Therefore,  equation  (A-l)  reduces  to 

A  P  * 

and  we  conclude  that  the  pressure  difference  across  the  blade  is  directly 

~  /  —  / 

proportional  to  the  swirl  schedule  VXu  We  will  call  V  the  swirl 

J  r7 

schedule  or  the  loading  distribution. 


Appendix  B  :  The  Bound  Vorticity 

In  this  appendix,  we  will  show  the  relationship  between  the  bound 
vorticity  (vortices  distributed  on  the  blade  camber  lines  to  model  their 
presence)  and  the  swirl  schedule  (or  the  gradient  of  the  pitch  average  velocity 
defined  in  chapter  2). 


I® 


®y 


The  flow  is  assumed  to  be  incompressible  and  inviscid,  and  the  far 
upstream  flow  is  assumed  to  be  uniform.  The  flow  is  thus  irrotational  and 
the  vorticity  must  lie  on  the  blade  camber  lines  (see  figure  above).  Therfore, 


we  require: 


JL  .  V*  = 

Moreover,  by  vector  identity 


(B-l) 


B-2 


To  satisfy  both  equations  (B-l)  and  (B-2)  and  the  conditon  that  the  vorticity 
direction  must  be  normal  to  the  x-y  plane  (2D  assumption),  we  can  write  the 
vorticity  field  as 

Jjl.  =  A  (*')  (  V*<  X  Vfr)  (B_3) 

where  VOO  is  the  "periodic  delta”  function  defined  in  Appendix  C. 

To  find  out  what  O"  is,  consider  the  circulation  around  path  C  shown 


in  figure,  we  may  write,  by  Stokes  theorem, 


/  i-*  -•// 

JC  JJA 

Note  that  the  line  integrals  along  path  (D  and  @  cancel  out  each  other  exactly, 
and  by  substituting  equation  (B-3)  into  the  right  hand  side  of  equation  (B-4), 
along  with  the  definition  of  the  pitch  average  velocity,  we  can  show  that 

I  ^  (X+AX,y)  -/  ^ 

f  +  rf+A 


(B-4) 


J  4«(-k) 


and  finally,  by  taking  the  limit  as  AX  -~>  0  ,  ^  in  equation  (B-3) 


defined  as 


Appendix  C  :  The  Periodic  Generalized  Functions 

The  "periodic  delta”  function,  the  "sawtooth”  function  and  the 
"smoothing”  functions  are  constructed  in  this  appendix. 

1.  The  "periodic  delta”  function 

The  "periodic  delta”  function  may  be  expressed  in  a  Fourier  series  of  the 

form 


h  =  -  to 


where  £  is  the  spacing  between  blade  camber  lines,  and  <<  represents 
the  blade  camber’s  surfaces. 

The  plot  of  SD  vs  o(  looks  like 
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2.  The  "sawtooth”  function 

The  ” sawtooth”  function  may  be  expressed  in  a  Fourier  series  of  the 

form 


Its  properties  are: 


-  It  has  first  derivative  related  to  the  "periodic  delta”  function  by 

$'(«)  =  J 

-  It  has  zero  average  between  fines 


(n  -i 


St*)  s  O 


-  It  has  a  jump  in  magnitude  of  A>  everytime  an  oC  surface  is 
crossed 


-  It  reduces  to  a  polynomial  form  when  (*-')  A  <  <<  < 

for  any  integer  rt  .  Its  polynomial  form  being 
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3.  The  ” smoothing”  functions 

The  "smoothing”  functions  are  defined  as 

Z  +\^oC 

where 


7 o  OxO 


/  ^  / 


Zirn, 


Their  properties  are: 


-  they  have  derivatives  of  the  forms 


H  =  (*)  V* 


-  they  have  zero  average  between  o<  lines 

n  A 

<=J*  «  O 

-  they  have  polynomial  forms  in  the  intervals  (rx— l)>d  <  K  < 

rv  A  for  any  integer  K  . 


L 


X|(*0  = 

2» 

- 

A* 

flC1 

It 

l 

2. 

s 

_  />*<< 
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a 

+-  ■ — 
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» 

K(«) 
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Appendix  D  :  The  "Physical”  Representations  of  /5  and  8 

In  this  appendix,  we  will  show  the  "physical”  representations  of 
Jk  and  S'  We  will  also  find  the  strength  of  the  vortices,  sources  and 
sinks  which  model  the  presence  of  the  blades. 

Consider  the  pitch  average  velocity  V  of  the  ’Finite-Thickness’ 

problem.  By  definition 


vc«)  =  j-f 


using  equation  (4.2-16)  for  Vx  ,  Y  can  be  shown  to  be 

J  *  (  *TX  V )  )  *3 

Since  VTx  and  s  are  the 
jk>  and  -  %  represent  the  x-component  and  y-component  of  the 
"imaginary”  flow  in  the  "blade”  region  respectively.  We  expect  these  variables 
to  be  proportional  to  the  blockage  distribution. 


(D-l) 


Ytu  are  the  gap  average  velocity  components, 


Consider  the  curl  of  V 

A* 

v  =  (v;,  -V)  id-2) 

In  chapter  2,  we  have  shown  that  if  the  vorticity  is  distributed  along  the 
blade  camber,  then  its  strength  is  related  to  the  gradient  of  the  y-component 
of  the  pitch  average  velocity.  Equation  (D-2)  is  indeed  the  case. 


Consider  the  divergence  of 

v.  1  =  (  V,  -  / )  P-3) 

-  '  / 

Since  the  flow  is  incompressible,  (  Vrx+  ft  )  represents  the  source/sink 
distribution.  Note  that  the  boundary  conditions  (  ■+*  jk  )  =  0 

(equation(4.2-17))  at  X  =  0  and  X  =  1  are  the  conditions  required  for 
the  blade  profile  to  close  there. 


Appendix  G  :  Numerical  Difficulties 

In  this  appendix,  we  will  attempt  to  understand  the  convergence  problem 
encountered  in  the  iteration  process  for  the  blade  camber  line  ~f-  when 
"partial  smoothing”  is  used.  We  study  the  mathematical  behavior  of 
and  %  . 

By  combining  AT,  &T  ,  £  and  ^  in  the  "partial  smoothing” 
forms  of  equations  (4.3-1)  and  (4.3-2),  we  arrive  at  the  following  boundary 
value  problems: 

%  ,  .'V  _  _  .  .  r'rH  .  „  / 


F(x)(f-i)X  +  X  +  X 

-  vT<  -  (iH  - 

with  boundary  conditions 

x'fo)  =.  o 
X'O)  -  o 

and 

F(*)(i'-i)r"+  f(xnYo-f'l  r'4  y 

*  -s  +  + 

with  boundary  conditions 


(E-l) 


(E-2) 


Y  5  -vT  +  f 
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X*)  = 


T(t) 

S£t) 


Since  -j-  itself  is  a  function  of  ^  and  ,  we  have  a  boundary  value 
problem  consisting  of  two  coupled  non-linear  second  order  ordinary  differential 
equations.  In  solving  them  iteratively,  we  expect  that  certain  difficulties  can 


The  strut  problem 


We  study  the  mathematical  behavior  of  and  ^  by  first  consider 

the  case  where  -j-  is  identically  zero  everywhere,  but  7~  is  finite.  This 
is  the  case  of  a  symmetric  blade  (or  a  strut).  In  this  case,  Y  is  also 
identically  zero  everywhere,  and  the  above  boundary  value  problem  reduces  to 


3CT)  ' 


3CT)  ' 

with  boundary  conditions 

A't6\  .  _  vY 


v' 

SO-)  Tx 


"  (-*) 
VDX  V 


(E-3) 


/$»  -  -  VT*  (o) 

/>(')  *  -  ^Tx  (0 


Consider  the  homogeneous  solution  of  the  above  differential  equation  by 
assuming  ^  to  be  of  the  form 

„  * 

A  =  e 

Substituting  it  into  equation  (E-3),  we  obtain 

-  :llt)  (if'.r)  +  i  =  o 

>c r> 


Assume  £* <  <  £■ ' 


Then 


Consider  the  case  where  £  <<  1.  Then,  by  definition 


3  (t)  a)  fr"  C 

t(t)  „  &  (-  tr ) 


and  thus 


It  can  be  shown  that  under  the  assumption  of 
that  e/e''  -  &(T)  and  the  assumption 


f A  <<  1,  we  can  say 

„  A- 

£  « c  b  justifiable,  except 


perhaps  near  X  =  0  or  1- 


Finally,  we  can  write 


As  an  example,  let 

T (%)  s  Z  A  (e>LOCK) 

where  BLOCK  b  defined  in  subsection  4.5.1  .  Then,  it  can  be  shown  that 


can 


With  the  assumption  that  <<  L  we  conclude  that 

be  a  highly  oscillating  function,  having  a  "natural  frequency”  of  the  order  of 


UOCK 


» 

) 


Finally,  we  look  at  the  amplitude  of  jl  By  definition,  i.e.  equation 


(4.2.13)  in  "partial  smoothing”  form 


t  ■  c'.-m')  -  i 

Therefore,  the  above  aaalysis  shows  that  jj  behaves  like 

/  &(*T) 

and  thus,  the  source/sink  distribution  (  Vjx  +  )  oscillates. 

A  computer  program  was  written  to  solve  the  differential  equation  (E-3) 
using  the  Chebyshev  collocation  technique  [7],  (8).  Numerical  results  show  that 

oscillates  with  a  natural  frequency  in  agreement  with  the  above  analysis. 

In  classical  aerodynamics  we  know  that,  for  a  smooth  blade  with  finite 
thickness,  the  source/sink  distribution  used  to  model  its  presence  should  exhibit 
like  a  "sine”  wave  with  a  singular  point  at  the  leading  edge.  Numerically,  in 
order  to  resolve  this  singular  point,  we  would  need  an  infinite  number  of  terms 
in  the  smoothing  series.  We  conclude  that  in  the  case  of  "partial  smoothing” 
(by  using  only  two  terms  in  the  smoothing  series),  we  are  unable  to  represent 
the  usual  source/sink  distribution.  However,  such  a  representation  is  also  not 
necessary. 

.f 

Accordingly,  we  will  call  (  +  A  )  the  "modified”  source/sink 

distribution,  and  show  that  it  can  be  used  to  produce  corresponding  blade 
shapes  satisfactorily.  Thus  we  9eek  a  practical  method  for  the  design  problem, 
using  "partial  smoothing”,  without  necessarily  having  to  go  into  more  extensive 
mathematical  development. 


The  loaded  blade  problem 


From  the  above  discussion  of  the  strut  problem,  we  expect  expect  both 
^  and  ^  to  have  oscillating  behaviors  when  the  blades  are  loaded. 
Two  iteration  schemes  for  -j-  were  investigated  using  "partial  smoothing”. 

We  note  that  equations  (E-l)  and  (E-2)  have  certain  symmetry.  Rewrite 


them  in  operator  forms: 


*<{*)  -  ■  S, 

*,{*]  +  Mx}  -  <** 


(E-4) 


where 


M  }  5  [F(*>  *']  II 

4  {) s  |y -  f('-  ]  1 1 

*  *  VT,  -  (Wi 

'  'I,*  (&l 

The  boundary  conditions  remain  the  same  as  in  equations  (E-l)  and  (E-2). 

Method  1 

In  this  method,  we  attempt  to  solve  equations  (E-4)  simultaneously  for 
j\ >  and  Y  using  the  Chebyshev  collocation  technique. 

We  express  ^  and  ^  as  Chebyshev  series  and  convert  equations 
(E-4)  into  matrix  forms: 


i*nn  -wh  = 


which  can  be  arranged  in  the  form 


M  | 

-M 

i 

/  \ 

W 

\  —  j 

k. 

1 - 

*7 

i — > 

M 

t  * 

f  —  % 

<  * 

X  and  I  are  solved  by  inverting  the  above  matrix  using  an  IMSL 
subroutine  [9]. 

Results  show  that  when  more  than  approximately  11  collocation  points 
are  used  in  the  calculation,  convergence  in  cannot  be  achieved.  and 

V  are  found  to  oscillate  and  their  Chebyshev  coefficients  fail  to  converge. 
When  around  51  points  or  more  collocation  points  are  used  in  the  calculation, 
the  iteration  process  diverges  rapidly.  Two  conclusions  can  be  made  from  the 
results  of  this  method: 


1.  f9  and  a  pocess  oscillating  behavior  as  predicted  by  the 

analysis  of  the  strut  problem.  By  using  more  than  51  collocation 
points,  the  numerical  calculation  tries  to  resolve  the  Gibbs 
phenomenon  at  the  leading  edge,  but  fails  to  do  so  because  of 
"partial  smoothing”. 

2.  The  iteration  process  can  diverge  rapidly  because  of  the  very 

nature  of  the  iteration  process.  We  note  that  even  though  we  are 
solving  X  and  Y  simultaneously,  we  are  in  effect  solving 

equation  (E-4)  iteratively  because  and  its  derivatives  are 
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updated  at  every  iteration.  We  can  therefore  look  at  the  iteration 
process  of  method  1  as  if  we  were  attempting  to  solve  for 

X  and  /  in  the  following  manner: 

a.  First  update 

-  *{Y-J  *  5," 

b.  Then  update 

K{r~!  -  -  Zlx-l's' 

where  is  the  iteration  level  in  the  iteration  process  for 

Since  operators  and  are  both  expected  to  have  normal 

mode  solutions  of  oscillating  behavior,  we  conclude  that  there  is  a 
chance  for  a  driven  reasonance  to  occur  during  the  iteration 

process  for  -f-  ,  which  can  lead  to  divergence  of  the  iteration 

scheme  itself.  We  think  that  this  is  indeed  the  case.  When  too 

many  terms  are  kept  in  the  Chebyshev  series,  higher  modes  are 
present  resulting  in  a  greater  chance  for  reasonance  to  occur. 

When  fewer  terms  (around  10)  are  kept  in  the  Chebyshev  series, 
we  are  staying  away  from  the  natural  frequency  of  the  operators 
and  resulting  in  a  stable  iteration  process. 

Method  2 

In  this  method,  we  use  the  iteration  process  described  in  section  4.3  . 
Derivatives  are  computed  numerically  using  two  methods:  Chebyshev 

collocation  method,  and  finite  difference  method.  Numerical  results  show  that 
similar  problems  as  in  method  1  were  encountered  here.  This  should  be 
expected,  as  observed  above,  because  iterations  are  being  used. 

We  decided  to  use  method  2  for  our  design  method  because  of  two 


reasons 


1.  method  2  is  much  more  efficient  than  method  1  (faster,  cheaper 
and  simpler).  In  method  1,  we  are  required  to  invert  a  matrix  at 
each  iteration  of  -£■  In  method  2,  we  are  required  to  compute 
derivatives  instead. 

2.  method  2  can  easily  be  modified  if  we  wish  to  keep  more  terms  in 
the  smoothing  series.  Equations  (E-l)  and  (E-2)  are  only  valid  when 
the  first  two  terms  in  the  smoothing  series  are  kept. 

Finally,  the  finite  difference  scheme  (central  difference)  is  used  to 
compute  derivatives  because  it  is  numerically  more  stable  than  the  Chebyshev 
collocation  method. 

A  computer  program  was  written  using  the  above  method.  It  is  found 
that  when  around  11  points  are  used  in  the  calculations,  convergence  in 
-f-  is  achieved  in  about  10  iterations.  When  more  than  around  20  points 
are  used  in  the  calculation,  fails  to  converge.  In  order  to  resolve  this 

problem,  we  propose  to  use  a  filter.  The  calculation  procedure  is: 


1.  when  more  than  11  points  are  used  in  the  calculation,  iterate  for 

using  a  "filter”. 

2.  when  11  points  or  less  are  used  in  the  calculation,  iterate  for 
'j-'  without  using  the  "filter”. 


Filtering  method 

Two  different  "filters”  are  developed  for  the  above  iteration  scheme: 


1.  -j~  is  filtered  using  a  least-squares  chord-wise  fitting  method 

[10],  The  combination  (<jy  +  T  ' — AT^  )  in  equation  (4.2-22)  is 
filtered  using  a  fourth  order  polynomial.  The  motivation  for  using 
a  polynomial  curve  fitting  method  is  that  we  expect  J-  itself  be 
represented  by  a  polynomial  of  low  order.  ' 


2.  the  pressure  coefficients  are  filtered  by  taking  the  average  of  the 
maximum  and  minimum  envelopes  of  the  C  _  curves.  The 


p  curves  are  filtered  only  near  the  leading  edge  region  (0  < 
X  <  -4).  The  envelopes  are  constructed  by  straight  lines  going 
through  the  maximum  and  minimum  points  of  the  Cp  curves. 
This  procedure  has  been  chosen  to  date  for  its  simplicity,  it  clearly 
admits  improvement  possibilities  near  the  leading  edge.  But,  till 
now  at  least,  this  approach  has  compared  adequately  with  known 
results  (see  Text). 


Appendix  F  :  Computer  Code  of  "Zero* Thickness”  problem 


c 
c 
c 
c 
c 
c 
c 
c 
c 

,c* 

'c 
c 
c 
c 
c 

REAL  X(lOl) .VMY(lOl)  .DVMY(lOl)  ,DDVUY(101).FM(101)  .*(101) 

1 .  DA  (101) ,  FNEV(lOl)  .F(lOl).DFOOl),  DDF  Cl  01)  .PLOT  (3. 101) 
COMMON /S. AO. A1 . BO . B1 . PI . I JX . X . UMAX . COR 
C 

C  READ  STATEMENT 

C 

READ(1 .  •)  S.  ALP1 .  ALP2, 1 JK.  ITER. EM 
WRITE (2. 50) S. ALP1 . ALP2. IJX 

SO  FORMAT  (5X, 'SPACING  S  *  ' . F6  3/5X. 'INLET  ANGLE  ALP1  *  ’ . F7  3/ 
15X. 'OUTLET  ANCLE  ALP2  =  \F7  3 /SX.  ’NUMBER  OF  POINTS  IJX  =  ' 
1.13///) 

C 

C  INITIALIZE  VARIABLES  FOR  CALCULATION  PURPOSE 
C 

NUAX=20 

PI=3  141592654 
RAD*S7  29577951 
TANl=TAN(  017453294«ALP1) 

TAN2=TAN(  0174S3294‘ALP2) 

XI JK=IJK-1 
DX=1  /XI JK 
XX=1  *DX 

sum=o 

DO  5  N=1.IJX 

XN*N 

SUM=SUM‘(1  /  (XN*XN)  ) 

5  CONTINUE 

C0R=1  6449341/SUN 
C 

C  COMPUTE  COMPUTATIONAL  LOCATIONS 
C 

DO  10  1=1. IJX 
J=IJX-1-I 
XX=XX-DX 
X(I)=XX 
PL0TC3. J)=X(I) 

10  CONTINUE 

C 

C  COMPUTE  INPUT  FOR  PARABOLIC  LOADING  CASE 
C 

DO  15  1  =  1. IJX 
i= I  JK  ♦  1  - 1 
XX=X(I) 

C0NST=6  • (TAN2-TAN1 ) 

VUT(I)=CONST‘(  S‘XX‘XX-XX‘XX‘XX/3  )‘TAN1 


PROGRAM  NAME  THIN  FOR 

MAIN  PROGRAM  FOR  INVERSE  DESIGN  OF  COMPRESSOR  BLADES 
2-D  INCOMPRESSIBLE.  INVISCID.  INFINITELY  THIN  THE 
LOADING  DISTRIBUTION  IS  OF  PARABOLIC  FORM 

INPUTS  S  -  SPACING 

ALP1  -  INLET  ANGLE  (DEGREE) 

ALP2  -  OUTLET  ANGLE  (DEGREE) 

IJX  -  NUMBER  OF  POINTS 
ITER  -  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED 
ERR  -  CONVERGENCE  CRITERIA  ‘ERROR* 
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DYlfir  ( I )  =C0NST«XX*  (1  -XX) 

ddvutu)=const*(i  -2  »xx) 

FU(I)=C0NST*(XX«»3/6  -XX*«4/12  )*TAN1*XX 
PL0T(1 .  J)=FM(I) 

DF(i)=ma) 

DDF(I)=DDVMY(I) 

D£N=1  ♦DF(I)*DF(I) 

F(I)-FM(I)«  08333333*S*S*DVMY(I)*(-1  •DF(I)»DF(I)/DEH) 

15  CONTINUE 

FNEf  (I  JK) =0 
F(I JK) =0 
C 

C  START  ITERATION  PROCESS  FOR  CAMBER  LINE 

,-C1 

0LDERR=100 
DO  1  KIT=1 . ITER 
C 

C  COMPUTE  A  AND  ITS  DERIVATIVE 
C 

DO  100  1=1.1 JX- 1 
DEN*1  ♦DF(I) *DF(I) 

A(I)=-DF(I)»DVMY(I)/DEN 

TERM1 = -DEN* (DDF ( I ) «DVMY ( I ) «DF ( I ) •DDVMTCI) ) 

TERM2=2  •DF(I)*DF(I) *DDF(I) •DVMY(I) 

DA(I)=(TERM1*TERU2)/DEN 

100  CONTINUE 
C 

C  COMPUTE  EDGE  VALUES 
C 

DEN0=1  ♦DFCI JK) **2 
DEN1=1  ♦DF(1)««2 

BO=DDVMT ( I JK) • ( 1  -2  *DF(I JK) •DFCI JK) /DEN0)/DE*0 
81=DDVUT(1).<1  -2  »DF ( 1 ) *DF ( 1 ) /DENI ) /DENI 
A0=-2  «DF  ( I  JK)  •DDVUY  ( I  JK)  •  C 1  -DF(I  JK)«DF(IJK)/DENO)/DENO 
At  =  -2  *DF ( 1 ) •DDVMY ( 1 )  •  ( 1  -DF(1)»DF<1)/DEN1)/DEN1 
C 

C  UPDATE  CAMBER  LINE  F(X) 

C 

XN0RM=-2  •FSUMCI JK. F) /S 

ERRUAX-0 

SUMERR-0 

XLOC=-l 

DO  101  1=1 . I JK 

DF(I)=V!TT(I)*  08333333»S*S«(-DDma)-A(I)*DDF(I)-DA(I)»DF(I)) 
1-2  •DFSUM(I.F.DF)/S 

101  CONTINUE 

DO  102  1=1 . I JK-1 

FNEA'(I)=FM(I )  ♦  08333333*S»S#(-DVUT(I)  -A(I)»DF(I))-2  •FSUM(I.F)/S 
1 -XNORM 

ERROR=ABS  ( FNEX  (I)-F(D) 

SUMERR=SUMERR ♦ ERROR 
IF (ERROR  GT  ERRMAX)XLOC=X(I) 

IF (ERROR  GT  ERRMAX) ERRMAX=ERROR 

102  CONTINUE 
C 

C  CHECX  FOR  CONVERGENCE  IN  F(X) 

C 

1RITEC2. 55) NIT  ERRMAX, XLOC 

55  FORMAT  (1  OX. ’ITER  »M2.'  . ERRMAX  =  ' .  F7  5 

1 ,  1  AT  X  =  '  .  F8  5) 
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AVC£RR=SUMERR/XI JK 

IFCERRUAX  CT  (1  UQLDERR)  )WRITE(2. S4)AVGERR 
S4  F0RMATC//2X.  ’•••••  ITERATION  SCHEME  DIVERGES  ""  •«•••’/ 
1  2X. '  AVERAGE  ERROR  *  \F9  S) 

IF (ERR VAX  GT  <1  1»0LDERR))C0  TO  998 
IF (ERR VAX  LE  ERR) GO  TO  998 
C 

C  UPDATE  VALUES  FOR  NEXT  ITERATION 
C 

DO  103  I=1.I«-1 
F(I)*FNE1(I) 

103  CONTINUE 

0LDERR=ERRMAX 

C 

1  CONTINUE 

C 

998  CONTINUE 

■'C  OUTPUT 
C 

SRITEC2.S0) 

60  F0RMAT(//T10. "X1 ,T2S, 'FN* ,T37. ‘FL01  ANGLE’ .TSS, ’F’ . T6? 

1. 'BLADE  ANGLE’/) 

DO  SOI  1*1. I JX 
J=I JK*1-I 
PLOT (2. I)=FNE1(J) 

ANGLFM=RAD«  ATAN ( VUT ( J) ) 

ANGLF =R AD* ATAN (DF ( J) ) 

•R I TE ( 2 . 6 l ) X ( J ) . FM (J) . AiCLFU . FNEI (J) . ANGLF 

61  FORMAT (S(SX. FI 0  S)) 

501  CONTINUE 

C 

C  CALL  JCF  PLOTTING  SUBROUTINE 
C 

CALL  QP ICTR (PLOT . 3 . 1 JX . QT < 1 . 2) . QX (3) . QLABEL (4) 
l.QTLABC  BLADE  [*2]  -  MEAN  STREAMLINE  [91)’) 

1 . QXLABC ’AXIAL  LOCATION  X’)) 

C 

STOP 

END 

C 

c  •• . •• . . . . 

c 

C  COMPUTE  FUNCTION  FSUM 
C 

FUNCTION  FSUM(I.F) 

REAL  F(lOl).XUOl) 

COMMON/S . AO . A1 . BO . B1 . PI . I JX . X . NMAX . COR 

REAL  LAKDM 

SUM-0 

DO  100  M-l.NUAX 

xv=m 

LAMDM=2  *P!*XM/S 

DELFO* (F(I)-F(IJK)) • LAKDM 

DELF1  =  (F(I)  -F(l) )  •LAL'DM 

TO=EXP  ( -LAMDM *X  (I))  /LAUDM**3 

T1*EXP(-LAMDM*<1  -X ( I ) ) ) /LAMDM* *3  . 

TERMO=BO*COS (DELFO) ♦AO*SIH(DELFO) 

TERU1 s-Bl *C0S(D£LF1 ) <A1 *SIN (DELF1) 

SUU*SUM*TO*TERMO*T1 •TEBM1 
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100  CONTINUE 

FSU‘J=  S<S<COR<SUU 
RETURN 
END 
C 

C  »»»..«« . . . ««» . »»««• 

c 

C  COMPUTE  FUNCT10R  DFSUN 
C 

FUNCTION  DFSUHO.F.DF) 

REAL  F(101).DF(101).X(101) 

COIQIOM/S.  AO.  Al.BO .  B1  ,PI .  I  JX.  X.  WAX .  COR 

REAL  LAMDM 

SUU=0 

DO  100  M=1.N!IAX 
XH~H 

LAUDU-2  <PI<XU/S 
OELFO- (F ( I ) -F ( I JX) ) • LAMDH 
DELFl<(F(I)-F(l))<LAlffiU 
COSO*COS<DELFO) 

COS1*COS(DELF1) 

SIXO-SIX(DELFO) 

SIM1=SIN(DELF1) 

EXP0=EXPC~LA1£DU<X(I) )  /LAUDU*<2 
EXP1*EXP(-LAHDM<(1  -X(I) ) ) /LAMDM«2 
TERMOO=-EXPO<  CB0<C0S0<A0<S1 NO) 
TERX01=DF(I)<EXPO<(-BO<SINO<AO<COSO) 
TERM 1 0= EXP 1< (-B1 <C0SI<A1<SI Ml) 

TERX1 1=DF (I )  <EXP1<  (B1  <SIN1<A1<C0SU 
SU'J=SUU*TEJtaOO»TERH01<TERM10<TERMll 
100  COMTIHUE 

DFSUU=  S<S<COR<SUU 

RETURM 

END 


C 


Appendix  G  :  Computer  Code  of  ”  Finite-Thickness”  problem 


. . . . . 

c  •  • 

C  •  PROCRAM  NAME  THICK  FOR  • 

C  •  MAIN  PROCRAM  FOR  INVERSE  DESICN  OF  COMPRESSOR  BLADES  • 

C  »  2-D  INCOMPRESSIBLE. INVISCID  NITH  FINITE  THICKNESS  • 

C  •  US INC  THE  CENTRAL  DIFFERENCE  METHOD  TO  COMPUTE  • 

C  •  DERIVATIVES  • 

C  •  • 

c  •••••••••••••••••••••••••••••••••••••••••••••••••••••*•••• 

C 

REAL  X(101) ,T(101) , DT(tOl) .DDTdOl) .  RS(lOl)  .RI  (101)  .  RJ(lOl) 

1  VMX(lOl)  .DVMX(lOl)  .DDVMXdOl) .  VMT(lOl)  .DVMY(lOl) 

1 . DDVUY (101) . FM(lOl) ,F(101) .FOLD(lOl) .DF(101) .DDF (101) 

1 .  BETA (101) .  DBETAOOI ) .  DDBETA(lOl) ,  DELT(lOl) .  DDELT(lOl) 

1 . DDDELT(lOl) . PHIHXB(IOI) . PHIHTB(lOl) .DPHIHX(tOl) .DPHIHT(lOl) 
l.XLOAD(lOl).XILOADdOl)  CPT(lOl) .  CPB(lOl) 

1 . TdOl) .  VXTdOl) .  VTT(lOl)  ,  VXB(lOl) ,  VTB(lOl) .  PLOT (6. 101) 

1  AAAI  (10,10).  OSCdOl).SMOSC(lOl).SMCFdOl).OPTdOl) 
l.A(lOl)  .DAdOl).DDAdOl).B(lOl)  .DB(101) 

REAL  LAKDM 

C0MM0N/C0MM1/S. AO, A1 . BO. 81 . PI 

C0XM0N/C0KM2/T .  DT .  DDT .  RS.  R  I.  R  J .  DVMX .  DWT .  DDVMX .  DDVMT .  FM 
READ ( S . * ) BLOCK . S . ALP 1 . ALP2 . I JX . UMAX . I TER . ERR , CA . CB . CC . CD . NOPT 
1 . IPOIER. BLADET . BLADES . PRES ST . PRESSB 

NRITEd.DBLOCK.S.ALPl. ALP2.UK.  ITER. ERR.NOPT.CA.CB.CC. CD 
1  FORMAT(///SX. 'INPUT  PARAMETERS ‘//TS. ‘MAX  BLOCKAGE  -  ’ . 

1F10  5/TS.  'SPACING  *  \F10  S/TS.  'INLET  ANGLE  =  '.F10  S/ 

ITS. -OUTLET  ANCLE  *  \F10  ill 
ITS  -NUMBER  OF  POINTS  IJK  =  ‘.13/ 

ITS. -MAX  NUMBER  OF  ITERATIONS  ALLOWED  *  '.13/ 

ITS. -MAX  ERROR  IN  FIX)  ALLOlED  ERRMAX  *  \F10  61 

ITS. -FILTERING  OPTION  *  • . 11//T7. -BLOCKACE  AND  LOADING  PARAMETERS1/ 

1T9.-A  *  ' .  F5  2/T9.-B  *  \FS  2/T9. 'C  *  \FS  2/T9.-D  *  \FS  111) 

C 

C  INITIALIZE  VARIABLES  FOR  COMPUTATION 

C 

PI *3  1415926S4 
RALPU  0174S3293«ALP1 
RALP2a  0174S3293*ALP2 
RAD*S7  29577951 
TAN1 *TAN (RALP1 ) 

TAN2*TAN (RALP2) 

VYN*  S»(TAN1*TAN2) 

V0NSET*(1  ♦VYN»»2)**  S 

AMPT-S»BLOCK/ (2  • ( (CA / (CA«CB) ) **CA) • ( (CB/ (CA’CB) ) '‘CB) ) 

XUK*UK-1 
DX*1  /XI JK 
a*  i  *dx 

SUM*0 

DO  666  N*1 . UMAX 
XN*N 

SUM=SUM*(1  / (XN«XN) ) 

666  CONTINUE 

C0R*1  644934068/SUM 
C 

C  COMPUTE  LOCATION 
C 

DO  20  1*1.  UK 

XX=XX-DX 

Xd)*XX 


20  CONTINUE 

IF (NOPT  E9  1)CALL  SHAI  (I  JK.  IP01ER  X.  AAAI) 

C 

C  CONFUTE  THICKNESS  INPUT 
C 

T< I JK) =0 
TC1)=0 

DO  21  1*2. UK-1 
XX*X(I) 

T ( I ) =AMPT • (Bt**CA) •  ( 1  -XX)**C8 

21  CONTINUE 

CALL  DERWUK.DX.T.DT) 

DT(1)=0 

CALL  DERIV(UK.DX.DT.DDT) 

C 

C  CONPUTE  INPUT  LOADINC 
C 

XL0AD(IJK)*O 
XLOAO ( 1 ) *0 
DO  22  1*2. IJK-1 
XX*X(I) 

XLOAD(I)*(XX**CC) • ( (1  -XX)**CD) 

22  CONTINUE 
C 

C  CENEBATE  T-CONPONENT  GAP-AVERACE  VELOCITT 
C 

XL0AD0*0 

CALL  XINTdJX.DX.  XLOADO.  XLOAD.XILOAD) 

ANPV* (TAN2-TAN1) /XILOAD ( 1 ) 

DO  23  1  =  1.  UK 

VNY(I)=ANPV*XIL0AD(I)»TAN1 

DVHT(I)=AKPV*XLOADU) 

23  CONTINUE 

CALL  DEBI V  ( UK .  DX .  DVNT .  DDVNT) 

C 

C  CONPUTE  FUNCTIONS 
C 

DO  2  1*1.  UK 
XX=X(l) 

KS(I)=  S*S-T(I) 

Bia)  =  -(S*S/12  )♦(  5*S*T(I))-(  S*T(I)**2) 

BJ(I)  =  -(S*S*T(I)/12  )•(  2S*S*TCI)**2)-(T(n**3/6  ) 

racn«  s«s/bs(i) 

DVNX ( I ) * (  5*S*DT(I))/(BSCI)«*2) 

DDVNX ( I )  =  (  S*S*DDTU)/BSa)**2)*(S»DTa)**2/HSCI)»»3) 
DF(i)=ma)/m(i) 

DDFa)*(m{I)*DVXT(I)-VNT(I)»Dm(n)/VUX(I).»2 
2  CONTINUE 

FN0*0 

CALL  XIHKUK.DX.FNO.DF.FN) 

C 

C  CALL  SUBROUTINE  TO  CONPUTE  EDGE  VALUES 
C 

DOTH* 2  •DDTIO/S 
RIO*RI  (1) 
mi*VNT(l) 

DDVNX1 * DDVNX (1) 

DDVNT 1 *ODVKT ( 1 ) 

FN1*FM(1) 

DF 1 *VNT ( I ) / VSX ( 1 ) 


F1*FV(1) 

CALL  EDGE  (ERR .  WAX .  S .  T  AR  1 .  T  AR2 .  DDTH .  R 1 0 ,  WT 1 .  DDVMX1 .  DDVMT 1 
1 . FX1 . FO. FI . DF0.DF1 . AO. A1 . BO . B1 . CCR) 

FQLDdJK)*FO 
F0LD(1)=F1 
DF(I JK)*DFO 
DF(1)*DF1 
Fd  JK)*F0 
F(1)*F1 
C 

C  CUESS  F  FOR  ITERATIOB  PURPOSE 
C 

DELF0--2  •FSUMCIJK. FOLD. T.RS.X.RMAX. COR) /S 

DO  130  1*2. I JK-1 

TERV1 =R J ( I ) »RS( I ) »RI ( I ) 

TERH2*(-Dmd)-DFd)*DVMX(I))/<l  •0F(I)*»2) 

DELF=-2  •  (TERMl«TERM2*FSUMd.  FOLD.  T.RS.X.RMAX.  COR)  )/S 
FOLDd)=FM(I)*DaF 
130  CORTIRUE 

CALL  DERIVd  JK.DX.  FOLD. DF) 

DF  d  JK) *DFO 
DF(1)*DF1 
C 

C  CALL  CUESS IRC  SUBROUTIRE  FOR  DBETA  ARD  DDELT 
C 

CALL  CUESS  ( I JX .  S .  X .  T .  DVMX .  DOVMX .  DBETA .  DDBETA .  DDaT .  DDDELT) 

C 

RRZTEC2. 10) 

10  FORMAT (SX.  ‘CUESS  IRPUTS'//TlO.  ’X1  .T2S.  *F’  .TIO.  'BETA*  ,T5S.  'DaT' 
1.T70.  -Da-  TBS.  ‘DTT '/) 

DO  11  JM. IJK 
1  =  IJXM-J 

DXX=Dra(I)*DBETA(I) 

DTY*-Dm(i)*DDani) 

XRITE12. 12)1(1)  FOLD(I)  .BETA(I)  DaTd)  .DXX.DTT 
12  FORMAT (6 (5X.F10  S)) 

11  CORTIRUE 
** 

C  ITERAT10R  PROCESS  FOR  CAMBER  LIRE 
C 

RRRR-0 

OLDERRMOO 

OLDERMOO 

DBETA(IJK)*-DVKXdJK) 

DBETA(1)*-DVMX(1) 

DDaT  dJK)*OVMT  (IJK) 

DDELT ( 1 ) *DVKT ( 1 ) 

DO  999  RRR*1 . ITER 
HHHR*MRflM*t 
C 

CALL  HOVOB (UMAX  I JX  FOLD . T. RS. X, PHIHXB . PHIHTB.COR) 

CALL  DERIVd  JK.DX.  PHIHXB.  DPHIHX) 

CALL  DERIVd  JK.DX, PHIHTB.DPHIHT) 

C 

DO  110  1*1. IJK 
DERM  ♦DF(J)**2 

Ad)*(DFd)*(’DVMY  (I)  »DDaTd) )  »DV)£X(1 )  •DBETA(I) ) /DER 
110  CORTIRUE 

CALL  DERIVd  JK.DX.  A,  DA) 

CALL  DERIVd  JK.DX.  DA. DDA) 


-112- 


c 

CO  lit  1*1, IJK 
DEV*1  *DF(1)*«2 

8(i)*(0Dma)-0DDara)*A<i)*DDF(i)*2  *da  d )  *df  d )  > /dem 

111  COVTIVUE 

CALL  DERIVCl JX.DX.8.DB) 

C 

CO  112  1*2. I  JK-t 

TERV*DT ( I ) • (RS (I ) *RI (I ) *R  J (I ) ) /RS ( I ) **2 

DBETA  ( I)*RJd)*  (DDA  (I )  -08  (1)  *0F  CI)-B(I)  *ODF  (I ) )  /AS  (I) 

1*TERV* (DAd)  *8(1)  *DF(I) ) -DPHIHXd) 
DDELTd)*-RJd)*DBd)/RS(I)-TEAH*Bd>*DPHIHTd> 

112  com  RUE 

CALL  DERIV ( I JK . OX . DBETA . DDBETA) 

CALL  DERIV ( I JK . OX . ODELT , OODELT) 

C 

C  COMPUTES  CAUSER  LIVE 
C 

DO  1SS  1*1. I JK 

OSC ( I ) * ( -DVMY ( I ) * ODELT ( I ) -OF ( I ) • (DVHX ( 1 ) •DBETA d)))/(l  ♦DF(I)**2) 
OPTd)*OSCd) 

1SS  COHTIKUE 

IFIVOPT  EQ  1 ) CALL  SllOOTHdJK.  IPOBER.X.  AAAI  .OSC. SMOSC) 

C 

ERRVAX=0 

SUV£RR=0 

XL0C=-1 

DO  ISO  1*2. I JK-1 
TERV1=RJ(I)  »RSd)  *RI  (I) 

IFCVOPT  EQ  1 )  OPT  ( I )  =  SMOSC  CD 

DELF=-2  • (TERSll *OPT ( I ) •FSUH ( 1 . FOLD . T . RS . X . WAX . COR) ) /S 
F ( I ) *FM ( I ) •DELF-DELFO 
ERROR* ABS (F ( I ) - FOLD ( I] ) 

SUV ERR = SUV ERR « ERROR 
IF (ERROR  CT  ERRVAX)XLOC*Xd) 

IF (ERROR  CT  ERRVAX) ERRVAX* ERROR 
160  COHTIKUE 

CALL  DERIVCI JK.DX. F.OF) 

DF(I JK)*OFO 
DF(1)*DF1 
C 

C  (RITE  COHVERCEVCE 
C 

IF(VOPT  EQ  DCO  TO  16S 
VRITEC1 . 1S2) VVVV . ERRVAX . XLOC 

162  FORUATtlSX.  'ITERATIOV  •  ’ . 13.  ' . '.'ERRVAX  «\F7  6 

1  •  AT  X  *  \F8  S) 

IF (ERR VAX  CT  (1  3«0LDERR) ) VRITEt 1 . 163) AVCERR 

163  F0RUAT(//2X. '•••••  ITERATIOV  SCHEME  DIVERGES  ••••■/ 

1  2X. ’  AVERACE  ERROR  *  '  F9  S) 

IFCERRVAX  CT  (1  3»OLDERR))CO  TO  99B 
IF (ERR VAX  LE  ERR) GO  TO  998 
CO  TO  169 
C 

165  ERRAVC*SUVERR/XI JK 
VRIT£(1 . 166)  MM ,  ERRAVC 

166  F0RVAT(15X. 'ITERATIOV  I'. 13. . 'ERRAVC  *'.FT  S) 

IF (ERRAVC  LE  ERR) CO  TO  99B 

C 

C  UPDATE  VALUES 


169 


DO  170  1*2. 1 JX-1 
FOLDd)*F(I) 

170  com  HUE 

OLDERR=ERRHAX 

OLDER*ERRAVC 

C 

999  CONTIHUE 

C 

998  CONTIHUE 

C 

C  (RITE 
C 

WRITE (1.60) 

60  FORMAT C//T10. 'X' .T24. 'LOAD (I) * .T40. *T(X) ' .TSS. ’FN(X) ‘ ,T70. ’F(X) ’/) 
DO  SOI  1*1. I JK 

J=I JK-l-I 
FL0T(1 . I)*X( J) 

PL0TC2, I)*F(J) 

PL0TC3. I)*F(J)»T( J) 

PL0T(4.I)*F(J)-T(J) 

PL0TC5, I)=BLADET 
PLOT (6. I)*BLADEB 
TERM=AXPV*XLOAD ( J) 

•RITECl . 61 ) X ( J) . TERM.TC J) .  FM(  J)  ,F(  J) 

61  F0RliAT(S(SX.F10  S)) 

SOI  CONTINUE 

CALL  QPICTR(PL0T.6.IJX.ar(2.3.4.S.6).QX(l).9UBEL(14) 
l.QTLABC BLADE  SHAPE')  .QXLA8C AXIAL  LOCATION  X')) 

C 

WRITE(2. 70) 

70  FORMATC//SX. '0UTPUT'//T10. 'X'.T2S. 'F'.T40. 'BETA'. TSS. 'DELT' 

1.T70.  'DXX'.TSS.  'DTT'/) 

DO  71  J=1 . 1  JK 
I=IJX«1-J 

DXX*DVKX ( I ) »DBETA ( I ) 

DTT**DVMT(I) ♦DDELT(I) 

IF(NOPT  EQ  0)CRITE(4. 7S)DBETA( J) .DDELT(J) 

7S  F0R1!AT(2(2X.  F10  S)) 

(RITEC2. 72)X(I) . F(I) . BETA(I) .DELT (I) .DXX.DTT 
72  F0RXAT(S(5X.F10  W) 

71  CONTINUE 
C 

C  CALCULATE  PRESSURE  COEFFICIENT 
C 

DO  800  1*1. IJX 
Yd)*F(I)*T(I> 

800  CONTINUE 

CALL  VEL(I JK. UMAX, COR. X.T.F.DF. DDF. RS.RI  RJ.PHIKXB.PHIHYB 

i  .  vux .  om .  Dora .  m .  d  vut .  ddvut  .  dbeta  .  ddbeta  .  ddelt  .  dddelt 

l.VXT.VYT) 

DO  801  1*1. IJX 
Y(I)*F(I) »S*T(I) 

801  CONTINUE 

CALL  VELdJX.NUAX. COR. X.T.F.DF. DDF. RS.RI  RJ. PHIHXB. PHIHTB 
1 .  VUX .  Dm .  DDVUX .  VUT .  DWT .  DDVUT .  DBETA .  DDBETA .  DDELT .  DDDELT 
1. VXB.VYB) 

WRITE (1 , 799) 

F0RUAK//T10.  'X*  ,T2S.  'Cp*'.T40.  'Cp-'  .TS3.  '  (d.lt»)P'/) 

DO  802  1*1. IJX 


799 


vT=(vxT(i)**2.mci)**2)**  s 
VB*(VX8CI)**2«VY8(I)**2)**  s 
CPT(I)*1  -CVT/V0BSET)**2 
,  *  CPB(I)*1  -(VB/V0RSET)**2 

802  COBTIBUE 

C 

810  DO  804  1*1 . XJX 

J*IJX«1-I 

DIFF=ABS(CPTCJ)-CPB(J)) 

PLOT (2. I)*CPT( J) 

PLOTO.  I)*CPB(J) 

PLOT  (4.  I)*PRESST 
PLOT(S.I)=PRESSB 

TRITEC3. *)X( J)  .CPTC J) .CPBCJ) .DIFF 
«R1TE( l . 80S) X ( J) . CPT ( J) . CPB ( J) . DIFF 
SOS  FORMAT (4 (SX.F10  S)) 

804  corrimiE 

CALL  QPlCTR(PL0T.S.lJX.QT(2.3.4.$).8X(l).QLA8aC14) 
l.QTLABC PRESSURE  C0EFFIC1EBT  Cp') . QXLABC 'AXIAL  LOCATIOR  X')) 
C 

C  (RITE  BLADE  SHAPE  TO  A  FILE  FOR  USE  IR  THE  RASA  DIRECT  METHOD 
C 

DO  $00  1=1 . I JX 
TRITECS.  *)X(I) 

900  CORTIMUE 

DO  901  1*2. IJX 
J=IJK*1-I 
TRITECS. *)X(J) 

901  CORTIMUE 

DO  902  1*1. I IX 
BSURF*F(I)-T(I) 

WRITECS. *)BSURF 

902  CORTIMUE 

DO  903  1*2. IIK 

J=IJX«1-I 

TSURF*F(J)«T(J) 

TRITEC6. *)TSURF 

903  CORTIMUE 
C 

STOP 

ERD 

C 

c  .«••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• 

c 

c  THIS  FUMCTIOM  COMPUTE  FSUMCI) 

C 

FUMCTIOM  FSUMCI. FOLD. T.RS.X.RKAX. COR) 

REAL  FOLD(lOl) .TC101) .RS(lOt) .XC101) 

C0UM0M/C0MM1 /S . AO . At . 80 , B1 . PI 

REAL  LAMDM 

SUM*0 

DO  100  M*1 . RMAX 
XU*M 

LAMDM *2  *PI*XM/S 

DELFO * C FOLD C I ) - FOLD C I IX) ) • LAMDM 

DELF1 *  CFOLD C I ) -FOLDC 1 ) ) • LAMDM 

TTT*SIM(LAMDM»T(I) ) /LAMDM'RSCI) *COS (LAMDM* TCI) ) 

T0*TTT*EXPC-LAMDM«XCI))/LAMDM**3 

T1*TTT»EXP(-LAMDM* (1  -XCI)))/LAMDM**3 

TERMO*BO*COS (DELFO). AO*SIM (DELFO) 


TERN1 *-Bl »COS (DELF! ) * A1 •SIN(DELFl) 
SU!i=SUM»I0«TEJU£O*Tl  •TERSI J 
100  CONTINUE 

FSUN*COR*SUH 

RETURN 

END 

C 

C  . . . 

,c» 

SUBROUTINE  DERIV(IJX.DX.F.OF) 

REAL  FUOl).DF(lOl) 

DF(IJK)*  S»(-F(Iflt-2).4  »F(1  JX-O-S  •F{I«))/M 
DFID*  S* (3  *F(l)-4  *F(2)»F(3))/DX 
DO  1  1=2.  UK-1 
DF ( I ) *  S»(F(I-1)-F(I«1))/M 
1  CONTINUE 

RETURN 
END 
C 

. . . . . . . 

c 

SUBROUTINE  XINT(IJK.DX.XINTFO.F.XINTF> 

REAL  F ( 101 ) . XINTF(lOl) 

SUV'XINTFO 
XINTFd  JK)=SUN 
DO  1  1*2. I JK 
J*IJK*1-I 

Sl»=SUN«  $.(F(J)«F<J.1))*DX 
XINTF(J)*SUN 
1  CONTINUE 

RETURN 
END 
C 


c 

C  SUBROUTINE  EDGE  EDGE  FOR 

C  THIS  SUBROUTINE  CALCULATE  EXACT  LE  AND  TE  VALUES 
C  IN  ORDER  TO  ACCELERATE  CONVERGENCE  IN  FIX) 

C 

SUBROUTINE  EDGE (ERR . NVAX  S.TAH1. TAN2 . DOTH .  RIO .  WT1 
1. CDViai, DDWT1. FN1.FO. FI. DFO. DFl.AO.Al.BO.il. COR) 
REAL  LAVDV 
PI  *3  141592654 
,  *  F0*0 

DF0=TAN1 

IC*0 

100  SUV0130 

suvu*o 

SUV02-0 

SUN22-0 

SUV03*0 

SUV33=0 

SUV00»0 

SUN«0 

DO  101  V*1.NVAX 
XV3V 

LAVDV<2  <PI<XV/S 
DELF3LAKDV*(F0-F1) 

CC3EXP(-LAMDM)/LAVDV«2 

SUVSIN3CC<SIN(DELF) 

SUVCOS3CC<COS  (DELF) 

SUV013SUV01<C0R*SUVSIN 
SUN 1 1 3SUV1 1 <COR  <SUVC0S 
SUV023SU1I02<C0R<SUVSIN/LAVDV 
SUU223SUU22*C0R*SUVC0S/LAVDV 
5UV033 SUV03<C0R<SUVSIN< LAVDV 
SUU33ZSUV33 <COR<SUVCQS< LAVDV 
SUX003SUVOO<C0R/LAVDV«3 
SUV 3 SUV* COR /LAVDV 

101  CONTINUE 
ZlzRIO*DDTH-(l  -DF1«2) 

ZZZz-(2  /RI0)<(SUH01m2-SUV11«2) 

Z2z(-  5<RI0<ZZZ) <DDTH 

Z22*(  S*RIO<ZZZ)*DDTH 

Z33 C4  /RIO) <SUV01<SUV11<DDTH 

SA1=((2  <DDTH<(SUV11<TAN1<SUV01))/RI0-DDVVX1) 

SA2-U-2  *DDTH* (SUN01  -TAN1*SUN11) ) /RI0-DDWY1) 

SA3*Z1*Z2 

SA4«2  <DF1<Z3 

SA223Z1-Z22 

DEN* SA3<SA22<SA4<SA4 

A1 3 (SA1 <SA22<SA2<SA41 /DEN 

Bl*  CSA2<SA3-SA1 <SA4J /DEN 

A0S2  <(-Al<SUNll-Bl<SUN01-l  ) /RIO 

B032  < (A1 *SUU01 -Bl <SUNl 1 -TAN1 ) /RIO 

XNORV3~BO*SUVOO* A1<SUV02<B1<SUV22 

FTE3FVI-BO<SUV22<AO<SUV02<B1<SUVOO-XNORV 

DPH I HT 3 -SUV <B0<SUV33<B1 -SUV03<A1 

DPHIHXsSUV<A0-SUV03*Bl-SUV33<Al 

RATIO3  (2  <DPHIHX-S<AO)  /  (DPHIHX«2-  25<S<S<A0) 

DFLE3  5<DPHIHT<RATI0 

XNUVTE3WT1<  75<B1<RI0<B0<SUV11<A0<SUV01 
DENTE31  <  75<A1<R10<B0<SUV01<A0<SUV11 
DFTE3XNUVTE/DENTE 


C  SUBROUTINE  NAME  CUESS  FOR 
C  THIS  SUBROUTINE  INITIALIZE  BETA  AND  DELTA 
C  TO  START  THE  ITERATION  PROCESS  OF  F(X) 

C 

SUBROUTINE  CUESS ( I JK . S . X . T . DVUX . DDVUX . DBETA . DDBETA , DDELT . DDDELT) 

REAL  DBETA(lOl)  DDBETA(lOl) . DOELT(lOl) .DDDELT (101) . DVMX(lOl) . DVMT(lOl) 
1 . DDVUX(lOl)  .DDWT(lOl)  ,X(101)  ,T(101) 

C 

PI=3  141S926M 
DO  1  1*1.  UX 
XX=X(I) 

TC0S-C0SC2  *PI*XX) 

TSIN=SIN(2  *PI»XX) 

DBETA(I)*-DVMX(I)»TCOS 

DDBETA(I)— 0DWX(1)«TC0S*2  •PI*DVMX(I)»TSII 
DDELT(I)*0 
DDDELT ( I) *0 
1  CONTINUE 


RETURN 


C  SUBROUTINE  NAVE  HOUQB  FOR 
C  THIS  SUBROUTINE  CALCULATES  THE  HOViJSENEOUS 
C  TER2S  PHIHXB  AND  PHIHTB 
C 

SUBROUTINE  HONOBOISIAX. I JK. F. T.RS . X. PHIHXB . PHIHTB. COR) 

REAL  F(101).X(101).T(101).RS(101).  PHIHXB  (101) .  PHIHTB  (101) 
COUVON/COUHl/S . AO . A1 . BO . B1 . PI 
REAL  LAVDU 
DO  200  1=1 . I JK 

suux*o 

SUUT=0 

DO  100  11=1. HUAI 
XV=N 

LAVDV=2  «P1»XM/S 
DELF0= (F ( I ) -F ( I JK) ) •LAMDN 
DELF1= (F(I) *F(1) ) *LAHDN 
XSIN=SIN(LAVDN*T (I) ) 

EXP0=EXP ( - LAVDN»X ( I ) ) 

EXP1=EXP(-LAKDN*(1  -X(I))) 

SINO=SIN(DELFO) 

SINUSIN(DELFl) 

COSO=COS(DELFO) 

COS 1 =COS (DELF 1 ) 

TO=XSIN*EXPO/LAVDH**3 

T1=XSIN«EXP1/LAVDN**3 

SUVX=SUia»TO»(BO»SINO-AO»COSO)*Tl«(Bl»SItl»Al«COSl) 

SUVT=SU!rfTO»(-BO«COSO-AO*SINO)*Tl»(-Bl»COSDAl»SINl) 

100  CONTINUE 

PHIHXB ( I ) =COR«SUVX/RS ( I ) 

PHIHTB ( I ) =COB.SUNT/RS { I ) 

200  CONTINUE 

RETURN 


C  SUBROUTINE  NAME  VEL  FOR 
C  THIS  SUBROUTINE  COMPUTES  THE  VELOCITY  FIELD 
C  IN  THE  CASCADE  RESION  (RECION  2) 

C 

SUBROUTINE  VELd  JX.  MUX .  COR ,  X.  T.  F.  DF.  DDF.  RS,  RI  .R J.  PHIHXB.  PHIHYB 
1 . VMX ,  Dm ,  DDVMX ,  VMT .  DVMY .  DDWfT .  D8ETA .  DDBETA .  DDELT .  DDDELT .  VX .  VY) 

REAL  X(101)  .TUOl)  ,F(101) . DFC101)  .DDF (101) . RS(lOl) . RI  ( 101) . RJ(IOI) 
l.PHIHXB(lOl)  .PHIHTB(101),m(101)  .DmClOl)  .  DDVMX(lOl)  . VHY(lOl) 
l. DVMT (101) .DDVMT (101) .DBETA(lOl). DDBETA(lOl). DDELT (101)  .DDDELT(lOl) 
l.VXdOl).VT(lOl).  TVX(lOl)  .  TVT(lOl) .  AX(101)  .BK101) 
l.PHIHXdOl)  .PHIHY(lOl) 

COMMON/ C0MH1/S . AO. A1 . BO. B1 . PI 
C 

CALL  HOMO (NM AX. I JK.F.X.Y.PHIHX.PHIHY.COR) 

C 

AX(IJX)=AO 
AX (1) SA1 
BXd  JX)*B0 
BX(1)=B1 

DO  200  1=2.1 JK-1 
DEN=(1  ♦DF(I)««2)»*3 
Xl  =  2  «DF(I)«(1  ♦DFd)»«2) 

X2=l  -DF(I)**4 
X3-DDF(I)»(1  -3  *DF(I)»»2) 

X4=DF ( I ) *DDF (I ) • (3  -DF(I)««2) 

DXX=DVUX(I ) ♦DBETA(I) 

DDXX=DDVMX (I ) ♦ DDBETA ( I ) 

DYT=-DVMY(I) ♦DDELT(I) 

DDYY* - DDVMT (I) *DDDELT ( I) 

AX (I) - (XI •DDYY*X2*DDXX»X3«DYY-X4«0XX) /DEN 
BX ( I ) = ( -X2»DDYY*X1 •DDXX»X4»DYT«X3*DXX)/DEN 
TVX(I)*( Drr-DFd)*DXX)/d  ♦0Fd)*»2) 

TVY ( I )  = (DXX^DF ( I ) *DYY) / ( 1  .DF(I)«»2) 

200  CONTINUE 
C 

DO  100  1*1. I JX 
ALPHA*! (I) -F(I) 

RST*  5*S- ALPHA 

RIT*-(S«S/12  )♦(  5«S«ALPHA)-(  5«ALPHA**2) 

CCC=RJd)/RS(I) 

TEP.UX=RSY*TVX(I)  ♦  (CCC'RIT)  *AX(I) 

TERMY=RST*TVY (I) ♦  (CCC^RIT)  •BXd) 

VXd )  =  VUX(I)  ♦PHIHX  (I)  * PHIHXB (I ) 'TERKX 
VT d ) = VUT ( I ) ‘PHIHY ( I ) -PHIHTB ( I ) ♦ TERMT 
100  CONTINUE 


C  SUBROUTINE  NAME  HOXO  FOR 
C  THIS  SUBROUTINE  CALCULATES  THE  HOVOGENEOUS 
C  TERNS  PHIHX  AND  PHIHT 
C 

SUBROUTINE  HOliO  (NHAX , I JK .  F.  X.  T . PHIHX .PHIHT . COR) 

REAL  F(101)  X(lOl) . T(lOl)  .PHIHX (101)  .PHIHT (101) 

COUNOH/COlffll/S ,  AO.  A1 .  BO .  81 .  PI 

REAL  LAKDM 

DO  200  1=1 . I JK 

SUHX*0 

suirr=o 

DO  100  Hzl.NVAX 
XN=H 

LANDN=2  *PI»XJI/S 
DELFO* (T (I ) -F ( I JK) ) •LAUDH 
DELF1=(T<I)-F(1))*LAMDM 
EXPO-EXP ( - LANDN*X  C I ) ) 

EXP1=EXP(-LANDN* (1  -X(I))> 

SINO-SIN(DELFO) 

SIN1=SIN(DELF1) 

COSO=COS (DELFO) 

C0S1*C0S(DELF1) 

T0=EXP0/LAHD¥«»2 

T1=EXP1/LAUDU««2 

SUUX-SUVX«T0< (-BO*SINO«AO*COSO) ♦T1»(B1*SIN1 »A1 >0081) 
suvT-suirr«To* (bo«coso«ao*sino) »ti» (B1 »C0S1 -A1 «SIN1 ) 
100  CONTINUE 

PHIHX(I)=COR»SUNX 
PHIHT(I)=COR«SUMT 
200  CONTINUE 


C  SUBROUTINE  MAKE  SMOOTH  FOR 
C  THIS  SUBROUTINE  SMOOTH  A  GIVEN  FUNCTION  F(X) 

C  USINC  A  MODIFIED  LEAST -SQUARE  SCHEME 
C 

SUBROUTINE  SMOOTH ( I JX. I ORDER. X.AI. F.SMF) 

REAL  XClOl).F(lOl). SKF(lOl) . AI <10. 10) .8(10) . SHCF(IO) 
C 

DO  400  J*l. IORDER-I 
SUM*0 

DO  500  1*1, I JK 

SUX*SUM*<F(I)-FCIJX>)*X(I)«J 
500  CONTINUE 

B(  J)-SUM 
400  CONTINUE 

8(I0RDER)*F(1)-F(I JX) 

C 

DO  100  1*1 . IORDER 
SUM-0 

DO  200  J=l. IORDER 
SUM*SUM*AI (I , J)*B(J) 

200  CONTINUE 

SMCF ( I ) =  SUM 
100  CONTINUE 

C 

DO  600  1=1. I JX 
SUM=F(IJK) 

DO  TOO  J*1 . IORDER 
SUM=SUM-SMCF( J) •(!(!) »»J) 

700  CONTINUE 

SMF  ( I )  =SUM 
600  CONTINUE 

C 


RETURN 


o  o 


:  SUBROUT t ME  MAKE  SMAI  FOR 
C  SUBROUTINE  USED  TO  CEMERATE  MATRIX  [  AI  ] 

USED  IM  THE  LEAST 'SQUARE  FITTING  FILTER 

SUBROUTINE  SMAI (I JK. IORDER.X. AI) 

REAL  AI  (10. 10)  .8(10)  ,«KAREA(20)  ,X(101)  .XM(H) 

C 

DO  109  I=2.2=I0ND£B 
SUM-0 

DO  200  J*l . I JK 

sus!*sua«xu)»»i 

200  CONTINUE 

xM(i*i)>sw 

100  CONTINUE 

C 

K*0 

DO  10  1=1 . IORDER 
00  20  1=1. IORDER 
Aid.  J)=XM(K*J) 

20  CONTINUE 

K=K*1 

10  CONTINUE 

DO  300  1=1. IORDER 
AI (IORDER. J)=l 
300  CONTINUE 

C 

Dl=-1 

CALL  LINV3F(A1 ,B, 1 . IORDER. 10. D1 .D2.IKAREA. IER) 


o  o  o 


c 

C  SUBROUTINE  NAVE  FILTER  FOR 
C  THIS  SUBROUTINE  FILTER  THE  PRESSURE 
C  COEFFICIENTS 
C 

SUBROUTINE  FILTERS!! JX.X.F.FFILT) 

REAL  X(lOl).FClOl)  .FT(lOl)  .XT(lOl) . FB(lOl) . XB(lOl) . FFILT(lOl) 
XIJK=IJX-1 
PI-3  141S926S4 
DX=1  /XI JX 
IJ=(IJK*l)/2 
C 

NB*0 
NT*0 
1=1 J-l 

100  1=1*1 

IFCI  CT  (IJK-1))C0  TO  999 
IF(F(I*1)  CT  F(I))CO  TO  100 
IFCI  ES  IJ)CO  TO  20 
NT=NT*1 
FT(NT)=F(I) 

XT(NT)=X(I) 

200  1*1*1 

IFCI  CT  <IJK-1))G0  TO  999 
20  IFCFC1«1)  LE  F ( I ) ) GO  TO  200 

N8*NB*1 
FBCNB)=FCI) 

XB(NB)=XCI) 

CO  TO  100 
999  CONTINUE 

FTCNT*1)*FCI JK) 

FB  CNB* 1 ) *FC I JK) 

XT  (NT*  1 )  *X  (I  JK) 

XB(NB*1)=XCI JX) 

C 

np=nt*i 

CALL  FITCI JK.DX.NP.X.XT.F.FT) 

NP=NB*1 

CALL  FITCIJX.DX.NP.X.XB.F.FB) 

DO  700  1=1. I JK 
FFILTCD*  5»CFBCI)*FTCI)) 

700  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  FITCI JX.DX.NP.X.XTB.F.FTB) 

REAL  XC101) . XTBC101) . FC101) , FTBC101) . FFIT(lOl) 
1.XFILTC101) .FILT1C101) .FILT2C101) .AAAI (10. 10) 

XNST=XTB(1)/DX 
NST*I JK-IIFIXCXNST) 

NST1=NST 
DO  10  1*2. NP 

X1C=(FTB(I)-FTB(I'1))/(XTB(I)-XTB(I-1)) 

XB=FTB  Cl) *XN»XTB ( I ) 

XJ=XTBCI)/DX 
JJ*I  JK-IIFIXCXJ) 


20 


DO  20  N=NST.JJ 
FFIT  01)  =XN*X(N)  »X8 
CO NT! HUE 
NST=JJ 
10  CONTINUE 
C 

DO  30  1=1 . NST1 
FTB(I)=FU) 

30  CONTINUE 

DO  31  I=«ST1»1.IJK 
FTB(I)=FFIT(I) 

31  CONTINUE 
C 

N=0 

DO  100  I=NST1. IJK-2 
N=N*1 

XFILT(N)=X(I) 

FILT1 CN) =FT8(I) 

100  CONTINUE 

CALL  SltAI (N. 4 , XFILT. AAAI) 

CALL  SMOOTH (N . 4 . XFILT . AAAI . FILT1 . FI LT2) 
N=0 

DO  101  I=NST1. IJK-2 
N=N«1 

FTB(I)=FILT2(N) 

101  CONTINUE 
C 

DO  200  I=NSTl-3.IJK-3 

FTB  ( I )  =  3333*  (FTB(I*1)*FTB(I) *FTB(I-1)) 

2 00  CONTINUE 

DO  201  I=NSTl-3 . I JK-S 

FTB ( I ) =  3333* (FTB(I»2)*FTB(I)*FTB( 1-2)) 

201  CONTINUE 
C 

RETURN 

END 

C 

c  . . . 


